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Continuous  wave  and  pulse  techniques  have  both  been  used  to  study 
deuteron  NMR  lineshape,  resonance  frequencies,  and  relaxation  times 
T  and  T„  in  oil-water  dispersions  as  a  function  of  temperature  and 
composition.  The  systems  studied  contain  hexadecane,  potassium  oleate, 
hexanol  (or  pentanol) ,  and  water.  The  water/oil  ratio  determines  the 
phases:  W/0  <  0.6  for  microemulsions ,  W/0  =1.0  for  cylindrical  rod 
liquid  crystal,  and  W/0  £1.2  for  lamellar  liquid  crystal.  A  model  in 
which  the  water  molecules  rapidly  exchange  between  a  "bound"  state  on 
the  surfactant  interface  and  a  "free"  state  in  the  bulk  liquid  is 
highly  successful  in  explaining  many  of  the  observations. 

For  the  microemulsions  the  relaxation  rate  is  proportional  to  the 
reciprocal  of  the  radius  of  the  water  droplet  spheres  dispersed  in  the 

oil.  The  correlation  time  for  both  bound  and  free  molecules  is  on  the 

-12 
order  of  10    sec,  which  is  fast  compared  with  the  reciprocal  of  the 

resonance  frequency(10   sec).  The  electric  quadrupole  interaction  is 

averaged  to  zero  by  the  rapid  interchange  of  water  between  the  two 

states  and  diffusion  of  "bound"  molecules  around  the  spherical  surface 
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of  the  droplets.  The  high  resolution  data  show   that  the  chemical  shift 
difference  at  23  C  between  the  two  OH  peaks  of  the  hexanol  and  H„0  is 
approximately  0.045  ppm  and  therefore  the  lifetime  of  the  OH  is  about 
250  ms.  The  chemical  exchange  rate  becomes  rapid  when  the  temperature 
is  raised  to  60  C  and  thus  explains  the  anomaly  in  T   occurring  at  57  C. 

The  cylindrical  rod  and  lamellar  liquid  crystal  phases  exhibit  a 
static  quadrupolar  splitting  of  several  hundred  Hz,  which  was  most 
readily  observed  as  a  CW  powder  pattern  spectrum.  The  observed 
quadrupole  splitting  is  approximately  ten  times  smaller  than  is  expected 
from  the  two  state  bound-free  model.  The  correlation  time  for  slowest 
tumbling  of  the  bound  water  at  room  temperature,  0.32  x  10   sec,  is 
only  slightly  larger  than  the  characteristic  time  of  experiment, 
approximately  0.15  x  10   sec.  This  explains  the  large  degree  of 
motional  averaging. 

The  pentanol  emulsion  has  two  coexisting  features  of  the  spin 
echoes,  one.  for  isotropic  microemulsion  and  the  other  for  a  liquid 
crystalline  structure  with  T  =  T„  -  0.72  ms .  The  portion  having  a 
liquid  crystal  characteristic  is  temperature  dependent  and  is  presumed 
to  provide  paths  for  charge  carriers  such  as  potassium  ions,  which 
enhance  the  electrical  conductivity. 
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CHAPTER  I   INTRODUCTION 

I . 1  History  and  Significance  of  the  Subject 

The  membrane  structure  has  been  a  major  interest  in  biology 
for  a  number  of  decades.  For  many  years,  biophysists  and  physical 
chemists  have  used  nuclear  resonance  techniques  to  study  the 
structure  and  mobility  of  the  "head"  group  in  lipid  membranes. 
This  is  basically  a  water-surf actant-lipid  interface.  The  properties 
of  the  lipid  and  surfactants  separately  are  understood  fairly  well 
from  extensive  studies  of  their  chemical  shift,  mass  spectroscopy,  UV 
and  IR  light  absorption,  etc.  But  upon  addition  of  water  the  system 
becomes  much  more  complicated.  The  principal  difficulty  encountered 
in  the  use  of  proton  nuclear  magnetic  resonance  lies  in  the  extreme 
difficulty  of  resolving  the  resultant  multicomponent  spectra  into 
assignable  resonances [1] , [2 ] ,  especially  when  these  individual 
resonances  occur  over  a  wide  frequency  range (typically  several  kHz), 
and  overlap  each  other.  Deuteron  magnetic  resonance (DMR)  was  first 
used  in  1971  to  study  the  mobility  of  the  hydrocarbon  chains,  in 
isotopically   enriched  molecules  by  Oldfield,  Chapman  and  Derbyshire. 
As  a  simple  model,  they  studied  the  modulation,  by  cholesterol,  of 
the  quadrupole  splitting  in  a  synthetic  di(perdeuterio)myristoyl-Lci- 
Lecithin,  as  a  function  of  temperature.  Subsequently,  Finer[3], 

Smith[4] ,  Seelig[5]  and  their  collaborators  have  performed  many  DMR 

2 
studies  on  the  hydration  of  Lecithin  and  of  the  dynamics  of  the  H 


dispersions  spin  labelled  lipids  occuring  in  membranes  and  model 
membranes.  The  kinetic  properties  of  the  hydratiar  water  in  such 
"liquid  crystals"  are  gradually  becoming  clear  to  us. 

The  structure  and  physical  properties  of  thermotropic  liquid 
crystals  and  lyotropic  systems  are  also  among  the  major  interests 
of  the  petroleum  and  pharmaceutical  industry.  They  are  intimately 
related  to  the  problems  of  oil  recovery,  and  pollution,  and  the 
mechanism  of  action  of  detergent,  med  cine,  etc. 

The  microemulsion  systems  and  liquid  crystals  studied  in  this 
dissertation  are  composed  of  hexanol  (pentanol) ,  hexadecane, 
potassium  oleate  and  deuterium  oxide  (DO).  To  describe  the  properties 
of  the  water  we  adopt  a  two  state  bound-free  water  model,  which 
assumes  a  layer  of  water  near  the  surfactants  having  prefered 
orientation  and  restricted  motion  while  the  remainings  of  water  have 
bulk-like  isotropy  in  motion.  The  validity  of  the  model  can  be  checked 
through  DiiR  with  both  continuous  wave  and  pulse  techniques.  The 
correlation  time  of  bound  water,  phase  transitions,  probability  of 
being  bound,  the  isotropic  rotational  process  of  free  water  molecules 
and  anisotropic  rotational  process  of  bound  water  molecules  are 
expected  to  be  deduced  from  the  experiment. 

A  brief  description  of  the  microemulsion  and  liquid  crystal 
phases  cf  the.  system  are  given  in  the  following  sections  of  this 
chapter.  An  equally  brief  basic  theory  of  nuclear  resonance  and  its 
application  to  this  particular  system  are  provided  in  chapter  II. 
Experimental  method  and  procedure  are  described  in  chapter  III. 
Chapter  IV  gives  results  and  conclusions. 


I . 2   The  Microem u Is ions 

1.2.1  The  Surfactant  and  Inter  facial  Tension 

Oil  and  water  can  be  dispersed  in  each  other  at  the  microscopic 
but  not  molecular  level  by  introducing  mixtures  of  surface-active 
compounds  as  emulsifiers.  The  mixtures  of  surface-active  compounds 
are  called  surfactant  and  co-surfactant  depending  upon  the  degree  of 
interfacial  tension  obtained.  The  systems  studied  in  this  dissertation 
are  composed  of  hexadecane(oil) ,  deuterated  water,  potassium  oleate 
(surfactant),  and  alcohol (co-surf actant) .  The  alcohol  used  is  hexanol/ 
or  pentanol,  but  we  put  emphasis  on  the  hexanol  system. 

The  potassium  oleate  creates  the  interfacial  film.  The  spontaneous 
development  of  negative  interfacial  tension  is  caused  by  penetration 
of  hexanol  through  this  film.  A  schematic  representation  is  provided 
in  Fig. (1) j  Fig. (2)  and  Fig. (3) [6]  on  the  next  three  pages  to  show 
the  hydrophilic  and  hydrophobic  duality  of  potassium  oleate  and 
hexanol. 

1.2.2  The  Formation  of  Microemulsions 

Microemulsions  are  optically  clear,  thermodynamically  stable 
oil-water  dispersions,  obtained  in  low  water/oil  ratios.  Phase 
diagrams  showing  the  regions  of  microemulsion  formation  have  been 
obtained  for  a  number  of  systems [7].  Studies  of  these  systems  by 
light  scattering [8] ,  low-angle  X-ray  scattering[9]  ,  viscosity [10]  , 
electron  microscopy [11] ,  ultra-centrif ugation[12] ,  nuclear  magnetic 
resonance[13] [14] ,  infrared  spectroscopy[13]  and  electrical  resistance 
measurements [14]  have  been  reported.  Water-in-oil  microemulsions 
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where  ft  is  the  interfacial  tension  reduction  due  to 
addition  of  hexanol 


Fig. (1)  A  schematic  representation  of  the  spontaneous  development  of 
negative  interfacial  tension  as  a  result  of  penetration  of  an 
interfacial  film  of  potassium  oleate  by  hexanol  at  the 
hexadecane /water  interface  [6] 
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Fig.  (2)   The  decrease  of  interfacial  tension  due  to  the  surfactant 
hexanol  [6] 
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Fig. (3)   The  decrease  of  interracial  tension  due  to  the  surfactant 
potassium  oleate  [6] 


have  been  shown  to  consist  of  water  droplets  from  less  than  100  A 

to  several  hundred  A  in  diameter  surrounded  by  an  interfacial  film  of 

surfactant  and  co-surfactant  as  in  Fig. (4). 

The  effect  of  hexanol  and  potassium  oleate  on  the  average  volume 
of  water  droplets  in  hexadecane  are  shown  in  Table(l)[15]  and  Table(2)[16. 


Table (1)   Effect  of  Hexanol  and  Potassium  Oleate  on  the  Average  Volume 
of  Water  Drops  in  Hexadecane  [15] 


Interfacial  Composition 
Hexadecane (oil)   Water 


Average  Drop-Volume 


No  additive 


No  additive 


Hexanol 


No  additive 


0.2  ml/ml  of  oil  No  additive 
0.4  ml/ml  of  oil  No  additive 
0.6  ml /ml  of  oil  No  additive 


No  additive 
No  additive 
No  additive 
No  additive 


Potassium  Oleate 
0.02  gm/ml  of  water 
0.06  gm/ml  of  water 
0. 1  gm/ml  of  water 


2.8x10  2ml 


6.0xl0~3ml 
5.0xl0~3ml 
5.0xl0~'3ml 


8x10  ml 
8xl0~4ml 
4.4x10  ml 


Hexanol 


Potassium  Oleate 


0.025  ml/ml  of  oil0.06  gm/ml  of  water 
0.05  ml/ml  of  oil  0.06  gm/ml  of  water 
0.075  ml/ml  of  oil  .06  gm/ml  of  water 


1x10  ml 
6.6xl0_5ml 
Volume  <10~  ml 


Table (2)   Effect  of  Water  Contact  on  Water  Droplet  Diameter  in  the 
Hexadecane  Systems  [16] 


Aggregate 


W/0  (*) 


Droplet  Diameter  A 


Inverted  micelle 
W/0  microemulsion 
W/0  microemulsion 
W/0  microemulsion 


0.1 
0.2 
0.4 
0.6 


32 
46 

56 
72 


(*)  1ml  hexadecane  was  presented  in  system  which  contained  0.2  gm 
potassium  oleate  and  0.4  ml  hexanol. 
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Fig. (4)      Schematic  representation  of    the    structure   of 
the  m icroemulsions 


I . 3   The  Liquid  Crystal  Phase 

When  the  volume  ratio  of  water  to  oil  increases,  the  microeraulsion 
breaks  down  and  a  liquid  crystal  forms.  For  the  composition  W/0  =  1.0, 
a  rod  like  structure  [Fig.  (5)]  appears.  Normally,  the  cylindrical 
structure  has  a  diameter  roughly  100  A  and  a  length  of  1000  A.  Water 
is  contained  within  the  surfactant  membrane.  The  membrane  stability 
has  been  discussed  by  Shah[17]  who  concludes  that  a  surfactant 
composition  of  1  hexanol  molecule  to  3  potassium  oleate  molecules  gives 
the  most  stable,  association.  It  is  speculated  that  potassium  oleate 
molecules  occupy  the  corners  and  hexanol  molecules  occupy  the  centers 

of  hexagons.  This  is  why  we  choose  the  system  with  potassium  oleate/ 

3 
hexadecane  =  0.2  gm/cm  and  volume  ratio  of  hexanol/hexadecane  =  0.4. 

As  the  water  concentration  is  further  increased,  a  lamellar 
structure  [Fig.  (6a)]  forms.  Additional  water  merely  increases  the 
gap  between  layers  of  surfactant  [Fig.  (6b)]. 

Figure (7) [18]  shows  the  change  in  the  viscosity  of  microemulsions 
and  liquid  crystals  (cylindrical  and  lamellar  structures)  upon 
increasing  the  water/hexadecane  ratio.  The  sharp  increase  of  viscosity 
after  disorientation,  such  as  violent  shaking,  is  due  to  the  mutual 
entanglement  of  the  liquid  crystals  [Fig.  (8)].  After  disorientation, 
the  cylindrical  structure  does  not  increase  viscosity  as  much  as 
lamellar  structure  because  of  easier  sliding  due  to  circular  cross 
section.  The  lamellar  structure  results  in  more  resistance  among  domains 
after  disorientation.  Thus,  the  fact  that  the  sample  of  W/0  =  1.0  has 
less  viscosity  than  the  sample  of  W/0  =  1.4  is  well  understood. 

The  gaps  between  liquid  crystals  of  both  cylindrical  and  lamellar 
structure  are  similar  to  slits  of  a  diffraction  grating  which  cause 
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light  (visible)  interference.  One  therefore  can  observe  colorful 
birefringence  through  polarized  plates. 

Oil  spheres  are  formed  for  water /oil  ratio  greater  than  2.0. 
In  this  case,  oil  is  surrounded  by  a  surfactant  membrane.  The  hydrogen 
bond  association  [Fig. (9)]  between  carboxyl  and  hydroxyl  groups  from 
adjacent  droplets  causes  an  abrupt  increase  in  viscosity  [Fig.  (8)]. 
This  kind  of  electrostatic  interaction  is  so  strong  that  the  viscosity 
is  much  larger  than  that  of  liquid  crystals.  This  region  is  not  of 
interest  in  the  present  dissertation  study. 
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Fig. (5)   Schematic  representation  of  the  structure  of  the 
birefringent  cylinders 
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Fig. (6)   Schematic  representation  of  the  structure  of  the 
birefringent  lamellae 
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Fig. (7)   The  changes  in  the  viscosity  of  microemulsions  and  liquid 
crystals  with  water/hexadecane  ratio  [18] 
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Fie  (8)   Disorientation  and  mutual  entanglement  results  in  higher 
viscosity  than  that  of  the  initially  well  oriented  liquid 


viscosity 
crystals . 
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Schematic  representation  of  the  oil  spheres  suspended  in  water 
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Fig. (9)   The  hydrogen  bond  association  between  carboxyl  and  hydroxy 1 
groups  from  adjacent  droplets  causing  an  abrupt  increase 
in  viscosity  for  water /hexadecane  ratio  of  2.0  to  3.5 


CHAPTER  II   BASIC  THEORY 

I I . 1  Nuclear  Magnetic  Resonance 

A  nucleus  with  spin  I  has  a  nuclear  magnetic  moment   |i=  Ynl, 
where  Y   is  the  nuclear  magnetogyric  ratio.  If  a  magnetic  field  B   is 
applied,  the  interaction  Hamiltonian  is 

}£=   -y  '  2  (II.  1.1) 

o 

The  nuclear  moment  precesses  about  the  applied  magnetic  field, 
and  the  quantum  mechanical  equation  of  motion  of  the  net  magnetization 
M  of  an  ensemble  of  spins  is  identical  with  the  classical  Newtonian 
law,  namely  the  change  of  angular  momentum  is  equal  to  the  torque 
applied.  Therefore 


~  =  yM  x  B  (II. 1.2) 

dt 


where  the  macroscopic  magnetization  M  is  defined  as  a  sum  over 
nuclear  magnetic  moments  of  all  nuclei. 

Since  the  nuclear  spin  I  is  space  quantized ,  the  allowed 
eigenstates  must  be  m  =  -I,  -1+1,  -1+2, . . . .1-1,  I.  In  thermal 
equilibrium,  the  distribution  of  population  on  each  energy  level 
is  governed  by  the  Boltzman  factor   exp  (-yfiB  /kT)  [see  Fig. (10)] 
to  guarantee  minimum  internal  energy. 
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£     TVN     and    Jl  =  exp(-YfiB  /kT)      (II. 1.3) 
m=-I  iJ   ,  o 

ra+1 

where  N   is  the  number  of  nuclei  in  the  eigenstate  designated  by 
magnetic  quantum  number  m,   N  is  the  total  number  of  nuclei  of  the 
ensemble,  and  T  is  the  temperature  of  the  lattice. 

The  statistical  description  above  strictly  requires  a  large 

20 
number  of  nuclear  moments  but  in  practice  N-  10   so  this  is  no 

problem.  Transition  between  states   is  induced  by  the  applied  rf 

field  with  equal,  probability  in  both  directions,  but  as  the  net 

population  of  the  lower  states  are  higher  than  those  of  the  higher 

energy  states  there  will  be  a  net  absorption  of  energy  from  the  rf 

source. 

A  steady  state  nuclear  resonance  experiment  can  be  interpreted 

as  the  result  of  two  competing  processes,  the  relaxation  which  tends 

to  establish  Boltzman  inequalities  after  nuclei  have  been  pumped  to 

higher  energy  levels  and  the  driving  electromagnetic  field  that  tends 

to  destroy  thermal  equilibrium  and  to  equalize  the  population  of  the 

different  energy  levels.  An  alternative  "classical"  description  is 

that  when  disturbed  from  equilibrium  the  nuclear  magnetic  moment 

process  about  the  magnetic  field  that  they  experience.  An  applied  rf 

field, orthogonal  to  the  main  applied  field,  will  interact  with  those 

precessing  nuclei.  An  oscillatory  field  can  be  resolved  into  two 

counter-rotating  components.  The  interaction  will  be  strongest  with  the 

component  rotating  in  the  same  direction  as  the  precessing  nuclei  and 

will  occur  when  the  precessional  frequency  is  the  same  as  the  rotation 

frequency  of  the  rf  field.  It  is  often  convenient  to  discuss  the 

situation  from  a  rotating  frame  of  reference. 
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II. 2  Relaxation 

The  macroscopic  magnetization  due  to  a  spin  system  in  a  large 
static  magnetic  field  relaxes  because  of  the  coupling  between  the 
microscopic  magnetic  moments  of  the  nuclear  spins  and  lattice  motions. 
As  the  molecules  move  about  in  space,  they,  of  course,  carry  the 
nuclear  spins  with  them  and  the  motions  of  the  nuclear  magnetic 
moments  generate  fluctuating  magnetic  field  which   is  experienced  by 
neighboring  spins;  this  time  varying  field  behaves  in  many  respects 
like  the  rf  field  generated  in  an  NMR  transmitter.  The  probability  of 
spontaneous  transition  is  negligibly  small  at  those  low  frequencies. 
Transitions  between  different  energy  levels  have, therefore, to  be 
induced.  They  may  be  induced  by  the  experiments  in  the  form  of  a  rf 
field  or  by  the  system  itself  by  a  random  motion  which  produces  an 
oscillatory  field  of  appropriate  frequency. 

Fluctuating  local  electric  fields  are  produced  also  by  molecular 
motion  and  where  the  nuclei  have  an  electric  quadrupole  moment,  this 
is  an  additional  usually  more  potent  source  of  relaxation.  To   describe 
the  return  from  non-equilibrium  to  equilibrium,  we  use  a  characteristic 
time  termed  T  ,  the  spin-lattice  relaxation.  T   is  often  described  as 
a  longitudinal  relaxation  time.  After  a  disturbance  the  nuclear 
magnetization  may  have  a  component  at  right  angles  to  the  applied 
static  field.  However  because  of  the  presence  of  various  local  fields, 
owing  to  the  fact  that  the  spins  are  actually  not  free  but  interact  with 
each  other  and  with  their  surroundings,  there  is  a  progressive 
dephasing  giving  rise  to  a  transverse  decay  of  the  magnetization.  We 
use  T„  to  define  this  "spin-spin  relaxation".  T?  is  also  called  transverse 
relaxation  time. 
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The  Bloch  equations  are  phenomenological  equations  developed  to 
describe  the  longitudinal  and  the  transverse  relaxation  mechanisms 
through  the  equation  of  motion (II. 1. 2)  of  the  nuclear  magnetization 
for  an  ensemble  of  free  spins.  In  a  laboratory  frame,  this  equation  is 


+  M  i  +  M  j    M  -M   /v 

§   =  Y  S  x  H  -   X   _,   y   -  -f-i  k  (II. 2.1) 

?         1 


where  i,j,k  are  the  unit  vectors  of  the  laboratory  frame  of  reference. 

CO 

Assume  that  the  applied  field  is  the  sum  of  a  d.c.  field  H  =H  =  -  — 

z   o    Y 

Wl 
and  of  an  rf  field  H\  of  amplitude   H  ,=  -  —  rotating  at  a  frequency 
1  1     y 

to  in  the  neighbourhood  of  to  .  This  field  will  be  the  more  effective 

o 

of  the  rotating  components  of  an  applied  field   H  =  2H.,  costot, 

x    1 

linearly  polarized  along  the  OX  axis  of  the  laboratory  frame.  In  the 
frame  rotating  around  H   at  the  frequency^  there  is  an  effective  static 

o  O 

field.  In  the  rotating  frame  the  observed  rate  of  nuclear  precession  equals 
yB  ±co   and  thus  the  effective  field  in  the  z  direction=  (yB  ±co)/y 
whereas  the  field  in  the.  i'  direction=  H ,  .  Neglect  of  the  counter- 
rotating  component  gives 


H   =  (H  +   -)kT  +  Hi'  (II. 2. 2) 

eff   o  y       1 


where  i',j',k'=k  are  the  unit  vectors  of  the  rotating  frame.  The 
equation  of  motion  of  (II. 2.1)  can  be  rewritten  as 


+               Mi'+Mi'  M  -  M   ~ 

dM    *•+       +        ,_     x     yr  z   o  ,  , 

dt   YU  X  eff;      T  Tx 

*      *  -> 

where  M  and  M  are  the  transverse  components  of  M  in  the  rotating 
x      y 


frame.  Therefore, 
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dM 
> 

dT 


-  +  (  03-.o)My 


dM. 
dt 


ox   r    l  z 


(II. 2. 3) 


dM       ^ 

■t —  =  w.M 
dt     1  y 


M  -  M 
z   o 


If  a  sufficiently  long  time  has  passed  to  allow  transient  effects 
to  decay  away,  we  can  obtain  a  steady  state  solution  by  setting 


dM    dM    dM 
x  _   y  _  2 

dt    dt    dt 


(II. 2. 4) 


Therefore, 


*    (  ««l0)  HI 
M  =  - 


2       2   2  2 
1  +  T„(  CO-OJ  )  +  Y  H.T -T, 


112 


H  T 
*      12 

M  =  — 


M  = 


1  +  T_(  oj-co  )  +  Y  H.T-T, 
2      o        112 

1  +  (  U-w  )2T^ 
o    2 

2       2   2  2 
1  +  T,(  W-U)  )  +Y  H-T^. 
2      o       112 


(II. 2. 5) 


We  can  reconvert  these  values  to  tne  laboratory  frame  by  using. 


M  +  iM  =  (M  +  1M  )  e 
x    y     x     y 


M  =  M  costot   -  M  sinojt 
xx  y 


(II. 2. 6; 
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M  =  M  sinidt  +  M  coswt  (I I. 2. 6) 

y    x         y 


The  components  of  the  magnetization  in  the  laboratory  frame  are 

functions  of  time  and  can  induce  in  a  coil  a  detectable  voltage  at 

the  frequency  w.  The  power  absorption  of  y)  =  -H«  -r-  only  occurs  in 

the  transverse  direction  and  has  its  maximum  value  at  u)=  id  by 

o 

i< 

contribution  of  the  term  M  . 

y 

Various  pulse  techniques  have  been  reported  in  the  literature [19] 

[20]  for  measuring  the  relaxation  times  T   and  T  .  In  this  work  the 

tt^t'vtt/2  pulse  sequence  was  used  to  measure  T1  and  the  tt/2'"^  t  ^  tt 

pulse  sequence  to  measure  T„ .  For  the  measurement  of  T   in  high  field, 

the  first  pulse  of  time  interval  At   ,    , tilts  down  the  magnetization 

pulse  1 

M  by  7T.  As  there  is  no  transverse  magnetization  only  longitudj.nal 

relaxation  exists.  After  time  T,  the  magnetization  vector  is  tilted 

back  to  the  x-y  plane  by  the  second  pulse  (  to  At   n    =  t7 / '2)  so  that 
J  J  o   pulse  I 

the  coil  can  pick  up  the  signal.  This  pulse  sequence  and  the 

measurement  of  T,  are  illustrated  in  Fig.  (11),  Fig. (12).  The  initial 

amplitude  of  the  free  induction  decay  signal  recorded  after  the  tt/2 

pulse  is  a  measure  of  the  net  component  of  the  magnetization  along 

the  z  axis.  In  Fig. (13),  the  first  tt/2  pulse  tilts  the  magnetization 

M  down  to  the  x-y  plane.  In  the  rotating  frame  coherent  with  the 
o 

Larmor  precession,  the  tilted  magnetization  spreads  due  to  the 

->- 
variations  in  H  over  the  sample.  The  second  pulse  applied  at  t  =  T 
o 

rotates  the  magnetization  by  tt  and  it  is  therefore  still  in  the  x-y 
plane.  The  components  of  the  magnetization  that  had  precessed  more 
rapidly  by  Ato=Y  AB  and  had  developed  a  phase  lead  of  A<f>=  Aid  T   now 
lag  the  mean  by  Acf> .  At  t  =  2t  ,  the  spread  magnetization  vectors 
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Fig. (13)   The  formation  of  an  echo.  Initially  the  net  magnetic  moment 
vector  is  in  its  equilibrium  position  (a)  parallel  tojzhe 
direction  of  the  strong  external  field.  The  rf  field  H.  is 
then  applied.  As  viewed  from  the  rotating  frame  of  reference 
the  net  magnetic  moment  appears  (b)  to  rotate  about  H  .  At 
the  end  of  a  tt/2  pulse  the  net  magnetic  moment  is  in  the 
equatorial  plane  (c) .  During  the  relatively  long  period  of 
time  following  the  removal  of  H. ,  the  incremental  moment 
vectors  begin  £o  fan  out  slowly  (d) .  This  is  caused  by  the 
variations  in  H   over  the  sample.  At  time  t=T  the  rf  field 
H   is  again  applied.  Again  the  moments  (e)  are  rotated  about 
the  direction  of  H  .  This  time  H   is  applied  just  long  etioug 
to  satisfy  the  pulse  condition.  This  implies  that  at  the  end 
of  the  pulse  all  the  incremental  vectors  begin  to  recluster 
slowly  (f).  Because  of  the  inverted  relative  positions 
following  the  TT  pulse  and  because  each  incremental  vector 
continues  to  precess  with  its  former  frequency,  the 
incremental  vectors  will  perfectly  recluster ed  (g)  at  t=  2~ ■ 
Thus  the  maximum  signal  is  induced  in  the  pickup  coil  at 
t=  2t.  This  maximum  signal,  or  echo,  then  begins  to  decay  as 
the  incremental  vectors  again  fan  out  (h) . [From  H.  Y.  Carr 
and  E.  M.  Purcell,  Physical  Review  _94_,  p630  (1954)] 
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overlap  again  and  an  "echo"  occurs.  The  envelope  of  echoes  generated 
by  a  series  of  tt  rf  pulses  is  governed  by  exp(-2t/T„)  [Fig.  (14)  ], 
that  is  if  diffusion  and  other  effects  are  ignored.  (   Note:  All  the 
pulse  intervals  are  assumed  much  shorter  than  T   and  T  .  ) 

II. 3  Deuteron  Relaxation  in  Water 

The  standard  Hamiltonian  representing  the  electric  quadrupole 
interactional]  is 

^Q  =  wfr-T>(  ^-I'Mnd^+if)]     (ii. 3.1) 

,  «A  2     2, 

where  eQ  =  <I,I \  *?   1  e.(  3  z.  -  R.)   I,I> 
1  i=l  i     i    i  ' 

The  two  parameters  used  to  define  the  electric  field  gradient  tensor 


eq  =  V        the  z  component  of  the  field  gradient 
zz 

V  -  V 

Ti  = "-    the  asymmetry  parameter  with   V   >  V   >   V 

i        y  j  J     r  i   z2 I _  I  yy1—  '  XX 

zz 


The  quadrupole  interaction  was  first  shown  by  Pound [22]  to 
dominate  the  spin  lattice  relaxation  in  water.  The  relaxation  times 
T..  and  T„  are  given  by  Bloembergen[2  3]  ,  calculated  by  standard 
time-dependent  perturbation  theory  applied  to  the  quadrupole 
interaction  term  (II. 3.1).  They  are 

_i  =  3  fLSQ   2  J^Q 8\ >  (  !  +  UL)    (ii. 3.2) 

T-    80^1   ;  {   ,    ,  2      2  +   .    ,     .      2   2  )    K       T   3  }         u     ; 
1  1  +  (0   T     1  +  A  W   T 
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i     7    2  n  9         5Tn  2xn  n2 

1_=  3_(e_aQ)2(      g q )(1  +  IL}   (TI.3.3) 

2  1  +  (0   T„    1  +  4  CO   T„ 

o  Q         o  Q 

Intensive  studies  have  been  undertaken  by  Eisenberg[24]  and 
Hindman[25]  on  the  quadrupole  relaxation  in  deuterated  water.  Even  in 
pure  DO  anisotropic  motion  does  exist  at  lower  temperature,  due  to 
cluster  of  hydrogen  bonded  molecules.  The  hydrogen  bond  influences  the 
activation  energy  and  in  turn  the  relaxation.  In  these  systems  the 
quadrupolar  relaxation  mechanism  is  generally  governed  by  two 
processes [25]  (I)  a  hydrogen  bond (dipole-dipole  pair)  breaking 
process  and  (II)  a  Brownian  isotropic  rotational  diffusion  process. 
These  two  processes  compete  with  each  other.  Process  (I)  dominates  at 
lower  temperature  and  process  (II)  dominates  at  higher  temperature; 
each  of  them  is  subject   to  a  different  activation  energy.  The 
relaxation  time  is  a  sum  of  two  Arrhenius  terms  as  shown  in  Fig. (15). 
For  simplicity,  we  only  focus  on  two  correlation  times  associated  with 
the  processes  mentioned.  If  exchange  of  DO  molecules  between  the 
hydrogen  bonded  lattice  and  a  isotropic  rotational  phase  is  rapid  a 
mean  correlational  frequency  will  result 


Tq-tJ  +  t*1  CH.3.4) 


and  in  turn  influence  the  observed  relaxation  rate, 


1=  1  ...  1  (II. 3. 4a) 

1    1 

T   and  T   are  the  T   relaxation  times  corresponding  to  the  correlation 

I    .   II 
times  T   and  T   . 


t 
(°c) 

I 

T«-12 

(   x  10        sec) 

5 

1.93 

20 

0.75 

40 

0.24 

60 

0.091 
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For  D„0,  the  averaged  quadrupole  coupling  constant  -  ^  over  a 
wide  temperature  range  is  about  258.6  kHz[25].  Correlation  times 
extracted  from  Hindman's  experimental  result  are  given  in  Table(3). 

Table  (3)   Temperature  dependence  of  various  correlation  time[25] 

11  t  -N   xa.   XI 

\17  TQ  (apparent)  =TQ+  TQ 

(  x  10  "sec)       (  x  10   sec) 

2.76  4.70 

2.03  2.79 

1.41  1.66 

1.03  1.12 

From  the  values  listed  in  Table(3),  we  obtain  two  correlation  times 

-12 
essentially  having  the  same  order  of  10   sec  in  all  cases  and 

therefore  we  can  omit  the  quadratic  term  in  the  denominators  of 

formulae  (II. 3. 2)  and  (II. 3. 3).  In  a  NMR  experiment  of  resonance 

frequency  6.5  x  10  Hz.  The  simplified  formula  will  become 

i--f(4»a)2Tq        whenn  =  0         (II. 3.5) 
■*■  >  * 

Activation  energies  due  to  two  different  processes  can  be 
calculated  directly  from  Fig. (15)  by  the  formula 


E   =  RT2  d  ln  Tl 
exp      — — (II. 3. 6)  ,R  is  the  ideal  gas  constant. 


For  the  convenience  of  comparison  with  experimental  results  on 
microemulsions ,  Hindman's  activation  energy  versus  temperature  graph 
is  included  in  Fig. (16).  The  figure  indicates  that  the  rotational  motion 
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Fig. (15)   The  contribution  due  to  (I)  hydrogen  bond  breaking 

process  and  (II)  isotropic  rotational  diffusion  process 
in  T  relaxation  of  pure  D  0[25] 
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Fig. (16)   Plot  of  apparent  activation  energies,  E&         (II. 2. 3.6) 
versus  temperature,  for  D„0  water [25] 
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for  the  D„0  molecules  relaxing  in  process  (II)  (high  temperature  path) 
is  isotropic  and  thus  rules  out  the  possibility  that  anisotropic 
rotation  is  a  factor  in  the  relaxation  process.  The  anisotropic  rotation 
will  be  reflected  in  differences  in  the  apparent  activation  energies 
for  relaxation  of  geometrically  non-equivalent  quadrupolar  nuclei  in 
the  same  molecule[26] , [27] . 

In  molecular  terms  water  is  normally  considered  as  a  mixture  of 
hydrogen  bonded  and  non-hydrogen  bonded  species.  The  nature  of  the 
mixture  may  be  continuous  or  discrete,  but  in  any  event  the  ratio  of 
non-hydrogen  bonded  to  hydrogen  bonded  will  increase  with  increasing 
temperature.  The  barriers  to  molecular  rotation  will  differ  between 
the  species,  and  hence  so  will  the  corresponding  activation  energies. 

I I . 4  Bound  Water  Model  for  Microemulsions 

Assume  there  are  two  types  of  water  [Fig. (17) ]  in  the  small  water 
droplets,  i.e.  (a)  free  water >  and  (b)  bound  water,  each  subject  to 
different  rates  of  molecular  reorientation.  The  bound  water  molecules 
have  their  prefered  orientation  determined  by  the  polar  groups  of  the 
surfactants  at  the  oil/water  interface.  If  the  radius  of  the  water 
droplet  is  r,  with  the  thickness  of  the  bound  water  layer   Ar,  the 
weighting  factor  of  bound  water  P,  is  given  by 

P  =  ^r2   Ar  m   3^r  (II. 4.1) 

b    4    3     r  u     ; 

■r  TTr 

Therefore,  the  weighting  factor  of  free  water  is 


3  Ar 
P  «  1  -  P  -  1  -  — -  (II. 4. 2) 

r        b         r 
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If  the  exchange  of  molecules  between,  bound  and  free  water  is  rapid, 

averaged  physical  observables  <  X  >  will  be  determined.  These  can  be 

the  relaxation  rate  —  ,  —  or  signal  height,  etc.  and  may  be  expressed 

1   2 
in  the  following  form 


<  X  >  =  Pu  X,  +  (1  -  P,  )  X,  (II. 4.  3) 

b   b         b    f 


and  for  the  microemulsions : 


<  X  >  =  Xf  +  3  Ar  (Xb  -  X  )  »  (II. 4. 4) 


If  the  observable  <  X  >  is  plotted  against  the  reciprocal  of 
radius,  i.e.  the  same  experiment  is  repeated  by  replacing  different 
microemulsion  samples,  a  constant  slope  is  expected  if  the  model  is  true. 
This    xf  can  be  obtained  from  the  intercept  of  the  line  with  the 
ordinate,  which  also  must  be  the  same  as  given  by  a  pure  D„0  data  at 
the  same  temperature. 

When  we  measure  the  height  of  the  free  induction  signal  as  the 

temperature  is  varied  an  abrupt  change  of  intensity  will  be  noted 

on  crossing  the  freezing  point.  The  loss  of  intensity  below  the 

freezing  point  is  caused  by  freezing  of  bulk  water.  The  intrinsic 

2 
quadrupole  constant   "   =258.6  kHz  causes  a  very  broad  line;  hence 

the  signal  amplitude  contributed  by  the  solid  DO  lattice  is  too 

small  to  be  detected  at  the  spectrometer  settings  used.  The  remainder 

of  the  signal  observed  is  totally  contributed  by  bound  water.  In  fact 

this  is  often  used  as  a  definition  of  bound  water.  The  "bound"  water 

is  that  water  whose  molecular  properties  are  modified  by  the  proximity 

of,  or  interactions  with  the  surface  to  such  an  extent  that  the 
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structure  of  the  water  is  not  competible  with  bulk  ice  and  this 
water  remains  unfrozen  below  the  bulk  freezing  point.  After  freezing 
the  relaxation  time  of  the  bulk  phase  is  reduced  and  exchange  with  the 
bound  phase  is  reduced  to  such  an  extent  that  exchange  between  bound 
water  and  bulk  ice  is  slow  or  the  new  characteristic  time  scale  and 
separate  signals  are  observable.  The  spectrometer  settings  can  be 
selected  so  that  only  one  phase,  i.e.  the  bound,  is  observed.  The 
indications  are  that  the  rate  of  molecular  motion  within  the  bound 
phase  is  approximately  the  geometic  mean  of  the  bulk  ice  and  bulk 
water  values.  Nevertheless,  the  bound  water  will  be  '"frozen" 
eventually,  when  the  temperature  is  further  reduced  and  the  rate  of 

molecular  motion  decreased.  Therefore,  the  change  of  intensity  at 

3  Ar 
the  freezing  point  gives  the  weighting  factor  P   of  — —  [  Fig.  (17)  ] 

at  that  temperature. 

Though  the  quadrupole  interaction  Hamiltonian  is  significant, 
we  do  not  expect  to  see  line  splitting  by  the  continuous  wave  method. 
The  bound  water  molecules  might  be  expected  to  have  a  prefered 
orientation  relative  to  the  droplet  surface.  However  the  occurance  of 
a  rapid  exchange  between  bound  and  free  water  and  the  small  size  of 
the  droplets  results  in  water  molecules  experiencing  regions  of  the 
bound  phase  at  all  orientation  relative  to  the  laboratory  >    I.e..    the 
magnetic  field  on  the.  NMR  timescale  thereby  averaging  out  the 
quadrupolar  spectral  features.  The  quadrupolar  interactions  are 
manifested  as  relaxation  mechanisms. 

The  exchange  of  nuclei  or  of  molecules  between  two  regions  has 
been  discussed  by  two  groups  of  workers.  Anderson  and  Fryer [28]  and 
Zimmerman  and  Brittin[29].  Anderson  and  Fryer  considered  the  transition 
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Fig.  (17)  The  bound  water  model  of  microemulsions 
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from  slow  to  fast  exchange  where  the  exchange  rate  was  compared  with 

the  molecular  reorientation  rates  in  the  two  phases,  whereas 

Zimmerman  and  Brittin  discussed  the  transition  when  the  exchange 

rate  was  compared  with  the  intrinsic  relaxation  rates  in  the  two  phases, 

At  temperature   to  the  high  temperature  side  of  the  1\  minimum  both 

the  spin  lattice  and  spin-spin  relaxation  rates  are  inversely 

proportional  to  the  molecular  reorientation  rate  v(=  —  )  as  in  the 

equations  (II. 3. 2),  (II. 3. 3)  and  (II. 3. 5)  although  the  proportionality 

relationship  is  more  general  than  for  relaxation  by  quadrupolar 

interactions.  It  should  be  noted  that  —  and  —  are  very  much  less 

1      2 
than  V  or  — .  If  the  nuclei  reside  in  the  separate  regions  or  phases 

TQ 
for  times  sufficiently  long  for  relaxation  to  occur,  separate 

relaxations  will  be  observed.  If  the  residence  lifetime  is  short 

compared  with  the  characteristic  relaxation  times,  population  weighted 

average  relaxation  rates  will  be  observed. 


P    P 

—  »/+-£  (II. 4.5) 

T    Tf   Tb 

Hence  a  small  fraction  of  bound  phase  having  a  greatly  increased 
relaxation  rate  can  dominate  the  observed  relaxation  rate  — . 


e.g.   if  T  =  1  s,   T  =  1  ms,   P  =  0.99  and   P  =  0.01 


1   0  99   0  01 

~  =  --—  +  — ~  =  0.99  +  10  and  T  a  100  ms 

10 


Zimmerman  and  Brittin  have  developed  expression  relating  the 

apparent  relaxation  rates  ~f     ,      — ,  and  P  '  and  P  '  to  the  exchange 

Jb     Tf 
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rate.  These  may  be  represented  diagramatically  as  follows 


(true  value)    T 


(true  value)    P 


SLOW 


FAST 


EXCHANGE  RATE 

Fig. (18)   The  relaxation  time  and  population  of  the  "bound"  and 
"free"  phases  influenced  by  exchange  rate 


In  developing  this  treatment  it  has  been  assumed  that  the  exchange 
process  does  not  itself  provide  a  relaxation  mechanism.  Because  spin- 
lattice  and  spin-spin  relaxation  times  can  have  different  values 
exchange  can  be  fast  for  one  measurement  and  slow  for  the  other.  It 
is  possible  to  generalize  the  treatment  to  apply  to  a  larger  number 
of  discrete  or  to  a  continuous  range  of  environments.  Therefore  fast 
exchange  in  the  Zimmerman  and  Brittin  sense  means  the  observation  of 
averaged  NMR  relaxation  rates  but  the  concept  of  separate  phases  has 
still  some  validity;  many  molecular  rotations  can  occur  in  each  phase 
before  exchange  occurs. 
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Anderson  and  Fryer  discuss  the  problem  where  the  residence 

lifetime  in  a  phase  becomes  comparable  to  and  ultimately  exceeds  the 

molecular  reorientation  rates  in  that  phase.  The  mathematical  treatments 

of  the  two  models  are  formally  similar  but  the  physical  consequences 

differ.  When  exchange  is  slow  in  both  models  nuclei  can  reorientate 

and  relax  in  their  separate  environments  before  exchange  occurs 

and  the  separate  environments  are  essentially  uncoupled.  As  —   ,  — 

1    2 
«  v   for  much  of  the  range  it  is  possible  for  exchange  to  be  slow 

in  the  Anderson  and  Fryer  and  fast  in  the  Zimmerman-Brittin  sense. 


LOG  REORIENTATION  &  RELAXATION  RATE 


LOG   TEMPERATURE 

Fig. (19)   The  relationship  between  reorientation  (or  relaxation) 
rate  and  temperature  (logarithmic  scale) 
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Measurements  of  reorientation  rates  e.g.  by  dielectric  techniques 
would  still  reveal  separate  phases  but  averaged  relaxation  rates  will 
be  observed.  As  these  are  inversely  proportional  to  the  reorientation 
rate  the  presence,  of  a  small  fraction  reorienting  slowly  will 
dominate  the  relaxation  as  discussed  earlier.  If  the  exchange  rate  is 
increased  to  become  fast  by  the  Anderson-Fryer  criterion  reorientation 
will  occur  at  a  rate  determined  by  the  population  weighted  average 
relaxation  rate  in  the  two  phases. 

V  =  P,  vc  +   P.  vu  (11.4.6) 

f   f    b   b 

-3 
Using  our  previous  values  P  =  0.99,  P  =  0.01,  V  =  10  Vf  ,  V 
°     t  a        b        b       r 

is  completely  dominated  by  the  abundant  rapidly  reorientating  a  phase 
and  hence  so  is  the  recorded  relaxation  rate.  This  again  assumes  that 
the  exchange  process  does  not  itself  provide  a  reorientation  mechanism. 

We  can  imagine  a  "thought  experiment"  where  we  increase,  the 
exchange  rate  between  two  phases  f  and  b  (for  convenience  both  may  be 
assumed  to  be  to  the  high  temperature  side  of  their  maximum  spin  lattice 
relaxation  rates),  without  at  the  same  time  modifying  their  intrinsic 
relaxation  and  reorientation  rates.  In  practice  this  is  impossible 
within  a  single  system.  An  increase  in  exchange  rate  is  normally 
accomplished  by  a  change  in  temperature  which  will  affect  the  intrinsic 
rates . 
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LOG  REORIENTATION 
&  RELAXATION  RATE 


LOG   TEMPERATURE 


We  are  envisaging  increasing  the  exchange  rate  along  the  line 


Fig. (20)   The  relationship  between  reorientation  (or  relaxation)  rate 

of  the  "bound"  and  "free"  phases  and  temperature  (logarithmic 
scale) 
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Fig. (21)   The  effect  of  exchange  rate  on  (A)  apparent  reorientation 
rates,  (B)  apparent  relaxation  rates,  (C)  apparent 
populations  of  NMR  (or  of  reorientation)  in  the  "bound" 
and  "free"  phases 
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The  core.  DO  has  only  isotropic  rotational  motion  with  rate  k  . 
The  bound  D„0  molecule  has  an  anisotropic  rotational  rate  k„.  We 
write  the  equilibrium 


k4 


loii 


DO  (bound)  +     '"  DO  (free)  .  (II. 4. 7) 

*"  k 

The  value  k   is  the  exchange  rate  from  free  deuteron  to  bound  deuteron. 

The  value  k.  is  the  rate  in  the  reverse  direction.  We  are  only  concerned 

2  2 
with  the  situation  CO  T  <<1,  which  is  true  from  the  argument  in  the  last 
o 

section.  Two  limiting  cases  to  be  considered  are  then:  (1)  the  case 

where  the  chemical  or  phase  exchange  is  fast  relative  to  the  rotation 

rate,  and  (2)  the  case  where  the  chemical  exchange  is  slow  relative 

to  the  rotation  rate  of  the  free  water  molecule,  i.e.,  the  relaxatic 

rate  is  limited  by  the  rate  of  breaking  polar  bonds. 

FAST  EXCHANGE:  k   and  k .  >>  k   and  k„ .  Tins  case  has  been  fully 
3       4     -t.       2 

treated  by  Anderson  and  Fryer [28]  with  the  result: 


^-  =  Q'[(l  +  K)/(k  K  +  k  )],  (II. 4. 8) 

1,2 


and  if  k?<<  k.K,  which  is  a  reasonable  assumption, 


¥~  =  Q'(  h+VT  }-  (II. 4.9) 

Tl,2      4  V 

where  K  =  k./k_  ,  and 
4   3 

Q'=  [(3/40) (21  +  3)/I2(2I  -  1)][(1  +  n2/3)(e2qQ/h2)]. 


We  note  that  the  assumption k  <<  k   is  based  on  the  premise  that 
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the  Debye-Stokes  equation  for  the  rotational  correlation  time, 


Tr-^#  (H.4.X0) 


where  T"]'  =  viscosity  and  a  is  a  molecular  radius,  should  be  obeyed  at 
least  approximately.  If  the  relaxation  were  controled  by  molecular 
motion  of  molecules  in  bound  sites,  i.e.,  k„>>  k  K, 


jp-  =  Q'(l/k  +  K/k  )  (II. 4. 11) 

1,2 

SLGW  EXCHANGE:  k,  or  k  >>  kn  and  k.  .  Anderson  and  Fryer  treated 
12     3      4 

tne  case  for  k  and  k  »   k   and  k  but  a  more  general  result  is  needed 
here  because,  as  noted  above,  k„  could  be  negligibly  small.  The  more 
general  T   equation  is  found  to  be 


jp-  =  Q'{[K/(l+K)](l+l/ki;  +  [l/(l+K)]fl/(k0+k/)] 

Tl,2  1  -      4 

+k,/(k  -k    )(k,+k    )}  (II. 4.12) 

3        12        I      4 


which   for   k  >>   k     reduces    to 


i-     =  Q'{[K/(l+K)](l/k   )   +    [l/(l+K)][l/(k  +k   )] 

1,2  L 

+k„/kn (k„+k.)}  (II.  4. 13) 

3      12      4 


The  subcases  are  then,  if  k„<<k. 

Z    4 


\r    =  Q'[(K/k  +l/k  )/(l  +  K)  +  1/Kk  ]  (II. 4. 14) 

1,2  H  L 


and  if  k  >>  k  and  k  and  K  not  «1, 
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±-  =  Q'  [(K/kx  +  l/k2)/(l  +  K)]  (11.4.15) 

1,2 


which  is  equivalent  to  the  Anderson  and  Fryer  result. 
For  K  «1 


±r    =   Q'/k  (II. 4. 16) 

1,2 


The  plotting  of  Fig. (22)  by  Anderson  and  Fryer  is  included  for 
future  data  analysis  of  the  significant  T„  change  of  the  microemulsive 
D„0  droplets  at  57  C.  This  graph  gives  the  change  of  relaxation  time 
due  to  chemical  exchange  of  fast,  intermediate  and  slow  passages  in 
the  presence  of  bound  water. 
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Fig.  (22) 


Change  of  relaxation  time  due  to  chemical  exchange 
of  fast,  intermediate  and  slow  passages  in  the 
presence  of  bound  water  [  Anderson  and  Fryer,  J.  Chera. 
Phys.  50,  3785  (1969)] 
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1 1 . 5  H igh  Field  Theory  of  Liquid  Crystal  Powder  Pattern 

Though  the  layers  of  the  liquid  crystal  are  ordered  within  a 
domain,  these  domains  are  distributed  randomly  throughout  the  sample. 
Thus  a  powder  pattern  [30] [31] [32]  is  obtained.  The  quadrupole 
interaction  therefore  must  be  weak  for  the  line  to  be  observable  in 
the  presence  of  this  broadening.  Only  first  order  splitting  need  be 
considered  in  the  present  case. 

&-^M+  ^Q  »here  >V<«m  (II'5'1) 

V  -  *  Ho  \ 

The  ~}-£      is  given  in  equation  (II.  3.1).  Since  the  rotation  of  the  DO 
Q 

-12 
molecule  is  fast  (about  10    sec),  the  electric  field  gradient  of  the 

OD  bond  is  averaged  around  the  molecular  axis  of  symmetry,  as  seen 

by  the  deuteron  nucleus.  For  the  isotropic  water  molecule,  the 

asymmetry  parameter  r\   is  zero,  but  is  not  for  an  anisotropic  water 

molecule.  The  anisotropic  DO  molecule  has  n  about  0.1  [25].  This 

small  asymmetry  parameter  can  be  ignored  in  high  field  theory  without 

causing  serious  error  because,  r|/2  appearing  in  the  quadrupole 

interaction  Hamiltonian  is  much  smaller  than  1.  The  nuclear  quadrupole 

interaction  Hamiltonian  of  (II. 3.1)  is  then  simplified  to  the  following 


KQ  =  A(3i;-I2)  (II. 5.2) 


where  A  is  a  constant  of  known  value. 

We  choose  xyz  as  the  strong  quantization  frame  (same  as  laboratory 
frame)  with  the  Oz  axis  parallel  to  the  static  field  Hq  and  XYZ  as  the 
principal  axis  frame  of  the  electric  field  gradient.  Since  the  system 


47 


is  cylindrically  symmetric,  we  can  choose  Oy,  OY  coincident.  Now, 

we  can  express  I„  in  terms  of  I  ,  I   and  0,  the  angle  between  the 
Z  x   z 

static  magnetic  field  H  and  the  principal  axis  of  the  electric  field 
gradient  of  the  molecule.  This  is  shown  in  Fig.  (23). 


xyz:  LABORATORY  FRAME 


XYZ:  PRINCIPAL  AXES 
OF  EFG 


Fig. (23)   Geometry  of  the  electric  field  gradient  principal  axes  at 

the  deuteron  nucleus  and  the  laboratory  frame  in  the  liquid 
crystal 


I„  -  I  sine  +  I  cosO 
Z    x        z 


Y    y 


(II. 5. 3) 


Also,  I  and  I  can  be  replaced  by  ladder  operators  I,  and  I 
x      y  + 


I  -  (  I  +  I  )/2 

x      +    - 


I  -  (  I  -  I  )/2i 
y       +     - 


(II. 5. 4) 
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Substitute  relation  (II. 5. 3)  and  (II. 5. 4)  into  (II. 5. 2). 

2 
}t Q(in  laboratory  frame)=  at?2I-1')  ^  2^3cOS  Q  "  1^3Iz  ~  I  "* 

+  |-  sinGcos9[I  (I  +  I  )  +  (I  +  I  )I 
2  z  +    -      +    -   ; 


3  2     2     2 

+  f-  sin  0  (if;  +  I  )] 

4  +    - 


(11.5.5) 


The  second  and  third  terms  do  not  contribute  matrix  elements  of 
the  diagonal.  The  first  order  perturbation  contributes  only  in  the 
case  of  m  =  0,  ±1,  therefore 


(o). 


YH 


E        =   -  Yhll  m  =   -  h   VTm  with  Larmor    frequency   VT    =  -z~r 
m  o  L  L        /.it 


E^=<m|KQ    |-> 


=  f^Yny     \   (3cos20   -l)(3m2-    1(1  +   1)) 


=  Tbvo(3cos20   -l)(m2-   1(1  +  l)/3) 
4       Q 


(II. 5. 6) 


where 


V  _  3  e  qQ 


Q   4  hi (I  -  1/2) 


For  the  deuteron,  I  =  1,  m  =  +1,  0 


E(l)=  e_aa  (3cos2q  _  1} 
1       o 


E(D=  _  e  aa(3cos20  -  1) 

o      4 


(II. 5. 7) 
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+  1 


-  _   V 


hV0  2 

w(0,-l)=  hv  +  — rH3cos  0-1) 

L         4 


hV0  2 

W(1,0)=  hV t^Ocos   0-1) 


(  KQ  -  o  ) 


m  Zeeman   splitting 


(  KQ  *  o  ) 


after  quadrupole  perturbation 


Fig. (2 4)  Energy  diagram  for  first  order  quadrupole  perturbation  of  the 
Zeeman  effect  of  the  deuteron  nucleus 


According  to  selection  rules  A  m  =  ±1,  the  transition  energies 
are  given  by 


W(1,0)  =  -  YfiH  +  -|  e2qQ(3cos23  -  1) 

O     o 


W(0,-1)  =  -  YfiH  +  1  e2qQ(3cos20  -  1) 
o    o 


(II. 5. 8) 


Calculation  of  the  spectral  shape  function,  which  is  a  product 
of  the  transition  probability  by  the  density  of  transition  at  each 
angle  0  ,  requires  the  transition  probabilities,  which  we  obtain 
directly  from  the  matrix  elements 

1/2 


<  m  1 1   m+1 
1  x1 


[(I  -  m)(I  +  m  +  1)] 


(II. 5. 9) 


<  m  1 1  I  m-1  >  =  -|  [(I  +  m)(I  -  m  +  D]1/2 
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1  =  1 


DASHED  LINE:  IDEALIZED  SHAPE 


SOLID  LIN 


3   2  A 


l-2x 


1-x 


1+x 


l+2x 


Fig. (25)   The  powder  pattern  resonance  line  for  I  =  1,  corresponding 
to  the  transitions   Am  =±1,  in  a  strong  magnetic  field, 
as  a  function  of  the  transition  energy.  The  ordinate  unit 
is  arbitrary.  The  nuclear  quadrupole  interaction  is  treated 
as  a  first-order  perturbation  of  the  NMR  transitions.  The 
figure  corresponds  to  a  nuclear  magnetic  resonance 
experiment  in  which  the  strong  magnetic  field  is  heJd 
constant,  while  the  oscillating  frequency  is  swept  through 
the  resonance. 
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therefore, 


P(1,0)    a    |<  l|    -2H   coswt  Yn  I    |    0>    |2   =   2y2fi2H   2cos2tot 

(II. 5. 10) 

P(0,-l)a  |<  0|  -2H  coscot  yfi  I  |  -1  >|   =  2y  H  H  cos  tot  . 


The  value  P  is  independent  of  6.  Thus,  the  shape  function  for  the 
transition  depends  on  9  only  through  the  density  of  transition  D  at 
angle  9.  For  an  isotropic  distribution  all  orientations  are  equally 
probable  and  the  number  density  of  nuclei  for  a  given  orientation 
varies  as 

d  cos9 


D(cos6) 


d  W 


(JJL.)-1 
M  cos9; 


2     W  ±  y^H       _  1 
£Xe  (— j-y-2- +  D  2  (II. 5. 11) 

a   e  qQ 


The  ±  signs  in  the  parenthesis  depend  on  which  one  of  W(0,-1)  and 
W(1,0)  is  chosen  for  W.  The  final  line  shape  is  determined  by 
superposing  the  shape  functions  of  both  transitions. 

The  idealized  shape  shown  in  Fig. (25)  by  the  dashed  line  has 
certain  striking  features.  For  values  of  W  corresponding  to  cos9  =  0, 
the  denominator  of  the  shape  function  vanishes  and  the  line  strength 
becomes  infinite.  It  is  also  to  be  noted  that  the  outsides  of  the  lines 
are  very  sharp,  while,  on  the  inside,  the  lines  merge.  Observed  spectra 
shown  by  solid  line,  may  be  expected  to  be  differ  from  this  idealized 
shape,  because  of  various  broadening  mechanisms,  such  as  lifetime  and 
dipolar  broadening.  In  addition,  the  effect  of  second  order  quadrupole 
perturbation  terms,  which  were  disregarded   in  our  discussion  of  the 
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line  shape-will  distort  the  shape  of  the  line  and  shift  the  positions 
of  the  minimum-without ,  however,  changing  the  distance  between  them. 

The  average  quadrupole  coupling  constant  of  the  deuteron  nucleus 
in  D90  has  the  magnitude  258.6  kHz  measured  by  Hindman[25]. 

The  simple  bound  water  model  of  lamellae  structure  is  shown  in 
figure (2 6).  If  the  water  layer  has  a  thickness  d=100  A  and  the  bound 
water  a  thickness  with  prefer ed  orientation  A=1.25  A,  the  observed 
quadrupole  constant  can  be  calculated  from  the  relation 


e\Q 


<  ^7-^  >  .   =  258.6  kHz  x  P  +  0  x  (1  -  P,  )     (II.  5. 12; 
h     obs.  b  b 


-2 


Only  bound  deuteron  water  contributes  to  the  quadrupole  splitting. 
The  weighting  factor  P.  for  bound  water  [33],  is  equal  to  A/d  =  1.25x10 

The  unbound  deuteron  water  with  fast  tumbling  frequency  (on  the  order 

12 

of  10  Hz,  is  much  faster  than  the  Larmor  frequency)  has  no  contribution 

due  to  quadrupole  splitting.  There  numbers  give 


2 

^-^  >  .   =3.23  kHz 
h    obs. 


(II. 5. 13) 


I 

A 

A  =1.25  X 

f 

c 

i  -loo  X 

OIL 
\\  \  \  \  \  \  \  \  \  \  \  \  \  \  \  .V.VAAAA. 

.'J  J  I  LL'J  LU  1 '-'  i  bi'J  i  - 
bound  WATER 

BULK  WATER 


OIL 

Fig. (26)   The  bound  water  model  of  lamellae  structure (bound  water 
has  prefered  orientation) 


CHAPTER  III   EXPERIMENTAL  METHOD  AND  PROCEDURE 

III.l   Sample  Preparation 

III. 1.1   Preparation  of  Chemicals 

The  chemicals  used  in  this  lesearch,  their  purities,  molecular 
weights,  and  the  commercial  source  are  listed  below: 

Table(4)   Components  used  for  sample  preparation 

Name  of  Component   %Purity   Molecular  Weight      Supplier 


Deuterium  Oxide  99.8 

(D20) 

Oleic  Acid  95 

Potassium  Hydroxide  85.9 

(KOH) 

1-Hexanol  99 


(C6Hl40) 

n-Hexadecane 

(C16H34} 
1-Pentanol 


99 


99 


20 


282 


56 


102 


226 


88 


Stohier  Isotope  Chemicals 


Fisher  Scientific  Co. 


Fisher  Scientific  Co, 


Fisher  Scientific  Co, 


Fisher  Scientific  Co, 


Fisher  Scientific  Co. 


Potassium  oleate  is  prepared  by  mixing  equimolar  amounts  of 
potassium  hydroxide  and  oleic  acid.  Since  solid  potassium  hydroxide 
does  not  readily  react  with  oleic  acid,  the  latter  is  mixed  with  the 
other  three  components (DO  ,  hexanol  and  hexadecane)  in  the  desired 
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ratio  first  and  then  with  potassium  hydroxide.  The  physical  properties 
of  the  dispersed  system  studied  are  highly  sensitive  to  change  of 
potassium  oleate. 

Since  potassium,  hydroxide  rapidly  absorbs  water  vapor  and  the 
amount  of  it  used  is  small,  it  was  necessary  to  quickly  and  precisely 

measure  the  weight  so  that  negligible  vapor  is  absorbed  during  weighing. 

-4 
The  resolution  of  the  balance  is  about  10   gram. 

The  samples  are  labeled  by  W/0,  the  volume  ratio  of  deuterium 

oxide  to  hexadecane  oil,  and  by  S/0,  the  ratio  of  surfactant 

3 
potassium  oleate  in  grams  to  hexadecane  oil  in  cm  .  All  the  samples 

have  a  fixed  ratio  of  the  cosurfactant  hexanol  or  pentanol  to  hexadecane 

oil  of  0.4  by  volume. 

The  prepared  samples  are  listed  in  the  following: 

Hexanol  system  (S/0=  0.2)  W/0  =  0.1,  0.2,0.3,  0.4,  0.5,  0.6,  0.7,  0.8 

1.0,  1.2,  1.4,  1.6,  1.8,  2.1 
Pentanol  system  (S/0=  0.2)  W/0=  0.1,  0.2,  0.4,  0.6,  0.8,  1.0,  1.2,  1.4 

Other  samples  were  prepared  with  changing  surfactant/oil  ratios  and 
a  fixed  W/0  ratio.  Unfortunately  most  of  them  were  not  stable  and  had 
a  separation  problem  during  trie  time  of  the  experiment.  These  samples 
were  all  rejected. 

III. 1.2   Sample  Sealing 

Since  the  sample  is  a  mixture  of  liquids,  small  amount  of 
dissolved  oxygen  molecules,  which  is  a  paramagnetic  impurity  having 

a  net  magnetic  moment  arising  from  unpaired  p  electrons  [34] [35] 

3 

roughly  10   times  larger  than  that  of  the  deuteron,  cause  a 

significant  contribution  to  the  relaxation  process.  In  order  to 
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exclude  the.  oxygen  molecules,  standard  freeze-pump-thaw  method  was 
employed.  Using  the  low  solubility  of  the  oxygen  molecules  at  cold 
temperature,  we  can  freeze  the  sample  with  liquid  nitrogen  first  and 
then  pump  to  evacuate.  The  next  step  is  to  thaw  the  sample  to  release 
the  remainder  of  the  oxygen  molecules.  Repeat  this  method  for  two 
cycles  and  then  seal  it.  This  procedure  represents  a  compromise. 
Prolonged  freeze -pump  thaw  cycles  may  well  cause  composition  changes. 

The  sample  container  is  a  Pyrex  tube  with  an  inner  diameter  of 
1  cm  and  an  outer  diameter  of  1.22  cm.  The  length  of  the  sample  tube 
is  about  3  cm.  Before  filling  with  sample,  the  tubes  were  cleaned  by 
aceytone  and  dried.  Sample  tubes  were  prepared  by  drawing  to  a  narrow 
necK  and  annealing  with  a  moderate  flame.  During  sealing,  caution  was 
exercised  to  heat  the  sample,  otherwise  the  chemicals  may  evaporate 
or  dissociate  even  though  most  of  the  sample  tube  was  submerged  in 
liquid  nitrogen. 

III. 2   Temperature  Measurement  and  Controlling  Apparatus 

III. 2.1  Cryostat 

The  cryostat  is  shown  in  Fig. (27).  After  locating  the  sample 
in  the  rf  coil,  the  inner  copper  and  intermediate  brass  chambers 
are  sealed  with  an  indium  wire  O-Ring  of  0.025"  diameter.  Helium 
exchange  gas  of  about  5  psi  is  admitted  to  both  chambers.  Though 
the  helium  exchange  gas  in  the  brass  chamber  was  latter  pumped  out, 
it  is  important  to  do  this  to  expell  water  vapor;  otherwise  the 
water  vapor  attached  on  electric  wire  may  cause  current  leaking 
especially  when  the  temperature  was  below  ice  point. 
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Fig. (27)  Cryostat  [ 
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For  experiments  below  room  temperature,  a  mixture  of  powdered 
dry  ice  (solid  CO  )  and  ethylalcohol  is  admitted  to  the  outside 
brass  cylinder.  Ethylalcohol  guarantees  good  temperature  uniformity 
of  the  heat  reservoir  at  -77  C.  Then,  the  helium  gas  between  the 

intermediate  brass  and  inner  copper  chamber  is  pumped  out  to  a  pressure 

-4 
of  10   torr  or  less  by  a  two-stage  vacuum  pump.  Even  though  the 

pressure  is  low  and  brass  is  a  poor  thermal  conductor,  slow  heat 

exchange  between  the  sample  and  environment  is  still  possible.  This 

slow  heat  exchange  facilitates  proper  cooling  of  the  sample.  The 

copper  chamber  has  excellent  thermal  conductivity,  thus  the  rate  of 

heat  exchange  between  heater  and  the  sample  is  much  faster  than  that 

between  environment  and  sample.  The  thermocouple  and  heater  are 

connected  to  the  feed  back  regulator  discussed  in  the  following 


II I. 2. 2  Thermometer 

A  copper-constantan  thermocouple  was  used  as  the  thermometer,  and 
is  located  at  the  top  of  inner  copper  chamber  [Fig . (27)  ].  A  mixture 
of  crushed  ice  and  water  was  used  as  reference  at  273.16  K.  Every  43  y.v 
thermovoltage  corresponds  to  one  degree  of  temperature  difference 
from  the  reference  point.  In  order  to  obtain  identical  temperature 
of  the  copper  chamber  and  the  sample,  measurements  snould  be  made 
during  heating  and  cooling  parts  of  any  thermal  cycle.  The  time  needed 
for  temperature  equilibrium  and  removal  of  the  temperature  gradient 
across  the  sample  varies  from  5  to  30  minutes  depending  upon  the 
temperature  difference  between  initial  and  the  final  states  and  the 
rate  of  heat  exchange  between  the  sample  and  the  environment. 


Ill . 2 . 3  Feedback  Regulator 

A  feedback  regulator  is  used  for  maintaining  the  sample  at 
constant  temperature.  Figure(28)  shows  the  block  diagram  of  the  feedback 
regulator.  The  three  amplifiers  in  the  system,  namely  the  FLUKE  823  AR 
DC-AC  differential  voltmeter,  the  DYMEC  2460  MI  amplifier  and  the 

Hewlett-Packard  467A  power  amplifier,  give  large  multifold  gain  of 

3 
approximately  10   so  that  even  a  small  difference  between  the 

temperatures  of  the  sample  and  control  value  is  sensitive  enough  to 

generate  a  considerable  quantity  of  additional  heat.  The  actual  and 

desired  thermal  potential  difference  between  the  sample  and  reference 

ice-water  mixture  may  be  expressed  as  kAt  and  is  compared  with  a 

voltage  E  .  The  value  V   generated  by  KEPCO  power  supply  (voltage 

regulated)  series  with  power  amplifier  IIP  467A.  The  difference  voltage 

is  amplified  to  produce  an  output  A(E  -  kAt).  It  is  necessary  to 

operate  the  amplifier  in  a  slightly  unbalanced  condition  so  that 

positive  and  negative  excursion  of  the  sample  temperature  from  the 

desired  value  can  be  compensated  by  decreased  or  increased  heater 

power.  If  operated  strictly  a  balance,  positive,  and  negative  voltage 

at  the  output  would  both  produce  heat  in  the  sample. 

The  unbalance  is  produced  by  the  insertion  of  a  "DC"  power  supply 

V  generated  by  the  KEPCO  in  series  with  the  balanced  HP  467A.  Even 
B 

so  if  the  output  of  the  amplifier  is  sufficiently  negative  the  system 
will  become  unstable,  a  high  temperature  will  generate  a  greater 
negative  voltage  (kAt  -  V  ),  and  a  greater  negative  current  will  flow 
through  the  heater.  Insertion  of  a  suitable  diode  prevents  the  negative 
current  flowing. 
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The  value  R^  is  the  resistance  of  the  sample  heater.  The 
distinguishable  temperature  differentiation  At  is  about  0.25  C. 


If  E  -  kAt  >  0  The  heater  current  turned  on  and  the  sample  heated  up, 
s 


E  -  kAt  <  0  The  heater  current  is  terminated  by  a  diode  between 
s 

the  HP  467A  and  KEPCO ,  and  the  sample  is  cooled  down. 

2 
The  required  value  of  V_,  is  such  that  V   /FL  is  larger  than  the 

cooling  rate  between  the  sample  and  environment  at  the  desired 

temperature.  The  regulating  procedure  are  listed  below 

(1)  Set  E  a  little  larger  than  kAt,  the  real  potential  developed 
in  the  thermocouple,  then  one  will  see  current  through  ammeter. 

(2)  Increase  the  gain  but  not  sufficiently  that  the  current 
exceeds  the  maximum  current  that  the  heater  can  tolerate. 

(3)  Adjust  VA  to  obtain  a  null  current  in  ammeter,  i.e.  the 
ammeter  is  balanced.  This  provides  proper  cut  down  current. 

(4)  Set  E   to  the  desired  value,  and  the  system  will  regulate. 
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I I I . 3  Nuclear  Res o nance  Measuring  Apparatus 

1 I I . 3 . 1   Circuit  Description 

Figure(29)  shows  a  block  diagram  of  the  probe  and  associated 
major  components.  The  output  impedance  of  the  transmitter  is  matched 
to  the  transmission  line  and  is  equal  to  the  input  impedance  of  the 
receiver.  The  impedance  at  resonance  of  the  sample  tank  circuit  is 
also  equal  to  that  of  the  transmission  line. 

The  operation  of  the  system  depends  upon  the  nonlinear  properties 
of  semiconductor  diodes.  During  the  application  of  a  high-powered  rf 
pulse  from  the  transmitter,  both  sets  of  crossed  diodes  in  Fig. (29) 
are  conducting  heavily  and  have  low  impedances.  Thus,  the  rf  pulse 
passes  through  the  diodes  at  A  without  appreciable  attenuation,  while 
the  diodes  at  B  prevent  the  rf  level  at  the  receiver  input  from 
rising  above  1  v.  Since  point  A  is  a  quarter  wavelength  from  B,  the 
effective  short  circuit  at  B  appears  as  an  open  circuit  at  A  and  thus 
does  not  affect  the  voltage  at  A. 

Following  the  rf  pulse  and  after  the  transients  have  decayed,  the 
diodes  stop  conducting  and  the  series  crossed  diodes  prevent  a  large 
impedance  at  point  A.  The  small  NMR  signal  voltage  (typically  about 
1  mv)  is  not  enough  to  lower  the  impedance  of  the  diodes,  and  thus 
in  the  absence  of  an  rf  pulse  the  transmitter  is  isolated  from  the 
system,  while  the  receiver  is  matched  to  the  sample  tank  circuit. 

It  should  be  observed  that  the  length  of  the  X/4  line  between 
points  A  and  B  is  not  specially  critical,  as  the  purpose  of  this  line 
is  to  allow  the  use  of  the  diode  short  at  B  without  presenting  a  low 
impedance  shunt  at  the  sample  coil.  Thus,  since  the  transformed 
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impedance  of  the  diode  short  need  only  be  large  compared  with  50  Q, 
within  the  bandwidth  of  the  tuned  sample  coil,  the  cable  need  only  be 
of  the  order  of  A/4,  and  the  use  of  the  quarter  wave  cable  does  not 
entail  any  significant  narrow-banding  of  the  system. 

A  more  detailed  block  diagram  is  given  in  Fig. (30).  The  rf  signal 
generator  (IIP  606B)  is  stabilized  by  a  HP  8708A  synchronizer.  High 
voxtage  pulses  are  provided  by  an  ARULAB  PG-6  50C  oscillator.  A 
programmed  pulse  generator  (Wang  Lab.  612AT)  controls  the  width  of 
pulses  and  the  delay  time  X      between  pulses.  The  resonance  frequency, 
pulse  width,  and  the  delay  time  can  be  checked  using  the  model  729C 
universal  time  counter.  Transient  signal   picked  up  from  the  sample 
coil  after  preamplif ication  (the  wide  band  amplifier  serves  as  power 
source  for  the  preamplifier)  is  sufficient  to  feed  the  phase  sensitive 
detector  (1-30  MHz  F+H  Instruments)  about  0.5  v.  However,  in  the  case 
of  spin-echo  measurement  of  T   the  stability  of  the  electromagnet  was 
not  adequate  to  use  phase  sensitive  detection,  and  diode  detection 
was  employed.  This  arrangement  is  shown  in  Fig. (31).  Signals  from 
the  phase  sensitive  detector  were  accumulated  in  a  Fabritek  1062 
signal  averager,  whence  the  resultant  could  be  plotted  on  an  X-Y 
recorder  (Honeywell  550). 

III. 3.2   Electromagnet 

This  is  a  Varian  V-3454,  9-inch  diameter,  Fieldial-regulated  unit, 
equipped  with  cooling  system.  The  drift  of  the  magnet  is  less  than 
240  parts  per  million  of  maximum  field  over  one  day  when  operated 
continuously.  In  other  words,  this  is  the  dominant  source  of  error  in 
the  maintainance  of  the  resonance  condition.  The  frequency  stability 
of  the  signal  generator  used  is  2  x  10   per  day. 
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Fig. (30)   Block  diagram  of  pulsed  NMR.  experiment  with  phase  sensitive 
detector 
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Fig.  (  31)   Block  diagram  of  pulsed  NMR  experiment  with  diode  detector. 
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11 1 . 3 . 3  Signal  Generator  and  Frequency  Synchronizer 

The  signal  generator  is  a  Hewlett-Packard  Model  606 B  unit, 
which  provides  a  signal  adjustable  from  0.1  microvolt  to  3  volts  and 
which  is  applied  to  a  50  ohm  resistive  load. 

The  frequency  synchronier  is  a  Hewlett-Packard  Model  8708A, 
which  is  designed  to  phase-lock  the  HP  606B  signal  generator. 

111. 3. 4  Pulse  Programmer  and  High  Power  Gated  Amplifier 

The  pulse  generator  is  a  Wang  Model  612  AT  transistorized 
programmed  pulse  generator.  This  unit  has  two  independent  channels, 
each  capable  of  delivering  twelve  pulses  which  are  collectively 
adjustable  in  duration  and  timing  as  follows:  the  pulses  in  the 
first  direct  channel  can  be  programmed  at  any  or  all  of  the  first 
twelve  equally  spaced  time  marks  established  by  an  internal  timing 
source;  the  pulses  in  the  second  delay  channel  can  likewise,  except 
that  in  addition  the  channel  as  a  whole  can  be  delayed  with  respect 
to  the  timing  source.  Since  the  period  of  the  timing  source  is 
adjustable,  a  flexible  pulse  program  is  available.  This  unit  has 
several  disadvantages:  The  change  of  delay  time  affects  the  pulsing 
rate.  There  are  only  twelve  delay  pulses  available  for  a  Carr-Purcell 
sequence.  Finally,  the  output  impedance  is  fairly  high,  and  the  output 
voltage  level  of  10  volts  maximum  is  marginal  for  the  requirements  of 
the  transmitter.  The  last  disadvantage  is  overcome  by  inserting  a 
high  power  rf  amplifier. 

The  high  power  rf  amplifier  is  a  Arenberg  Ultrasonics  Laboratories 
Instrument  Company  Model  PG  650-C  pulsed  oscillator.  The  unit  receives 
the  rf  signal  from  the  rf  signal  generator,  together  with  gate 
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rectangular  pulses  from  the  pulser.  The  transmitter  delivers  a  maximum 
rf  amplitude  of  about  300  volts  into  50  ohms. 

111. 3.5  Sample  Coil  and  Probe 

The  sample  coil  consists  of  16  turns  of  No. 22  double  cotton 
copper  wire  with  a  total  inductance  of  17yh. 

The  probe  (including  sample  coil)  serves  as  part  impedance 
matching  network  to  the  transmitter  and  receiver.  The  capacitors  of 
the  probe,  C.  ,  C„,  and  C,  are  tuned  so  as  to  maximize  the  observed 
signal.  The  circuit  diagram  is  given  in  Fig. (32). 

111. 3. 6  Receiver 

The  preamplifier  used  is  an  ARULAB  Model  PA-62  0-L.  This  is  a 
tuned  amplifier,  which  provides  voltage  gain  of  approximately  60  db 
with  a  low  noise  figure.  The  output  and  input  stages  each  has  nine- 
tuning  coil  caps.  These  are  used  to  balance  the  capacitance  existing 
between  grid  and  plate  of  the  vacuum  tubes  and  thus  limit  the 
recovery  time  and  optimize  the  system  [Fig. (33)].  The  frequency 
selectors  must  obviously  be  set  in  the  proper  range  for  sample 
resonance.  The  bandwidth  selector  can  be  used  to  adjust  the  resistance 
between  grid  and  ground  so  that  high  gain  with  narrow  bandwidth  or  low 
gain  with  broad  width  can  be  obtained  depending  on  the  frequency 
distribution  of  the  sample  signal.  The  well  tuned  narrowest  bandwidth 
is  about  200  kHz,  which  is  wide  enough  for  the  quadrupole  perturbation 
frequency  range  found  in  these  studies. 
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An  ARULAB  Model  WA  600  E  Wide  Band  .Amplifier  follows  the 
preamplifier  PA-620-L.  This  unit  has  an  rf  gain  of  about  60-65  db , 
flat  to  within  +3db  from  2  to  65  MHz. 

The  dead  time  is  about  1.5  x  10    sec.  The  gain  of  the  preamplifier 
is  about  180  and  the  total  rf  gain  of  the  receiver  is  about  10  . 

III. 3. 7   Phase  Sensitive  Detector 

The  following  is  a  qualitative  description  of  the  ring  diode 

phase  sensitive  detector.  For  a  more  rigorous  treatment  reference 

should  be  made  to  Arundell,L. , Clarendon  Laboratory,  Oxford  memo. 

The  assumption  is  made  that  V_  >>  V  and  that  the  diodes  when 

K      S 

forward  biassed  have  a  resistance  R  and  when  reverse  biassed  have  an 

infinite  impedance. 

Following  the  nomenclature  of  Fig. (34)  in  the  half  cycle  when 

V  >  V   the  circuit  becomes 
B    D 


Vg=  Vgsin(tat    4tj>) 


*L 


Applying  Kirchoff's  Laws, 
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2  vR  =  IXR  +  (Ix  +  II)  R  -  2  IXR  +  ILR 


(III. 3. 7.1) 


VR  "  Vs  =  %  +  (I1  +V  R  =  V  +  2  V 


(III. 3. 7. 2) 


2  vc 

•'•  -  2  v  =  3  I.R    or    I.R  = =-^ 

s      L  L       3 


(1II.3.7.3) 


In  the  other  half  cycle 


vD   =  V„sinijJt 


pnrTrmn 


v  =   V  sin(oJt  +   4) 
s        s 


As  before 


2  vR  =   IXR  +    (Ix  +   IL)    R  =   2    I  R  +   ILR 


(III. 3. 7. 4) 


VR  "  Vs   =    \ 


R     +    (I     +  I    )    R  =    I   R  +   2   I   R  (III. 3. 7. 5) 


2   vc 

and  -  2  v     =   3   IR  or  ITR  = ~ 

s  L  L  3 


(III. 3- 7. 6) 
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v    :  solid  line 
K. 

v  :  dashed  line 


Output 


h\ 


i 
I 
h-/ 

i  / 
1/ 


Ov 


If  the  assumption  is  made  that  the  switching  is  done  by  v   the 
putput  is  proportional  to 


-  V   ^   sin(0  +  d>)  d9 

TT   S     O  r 


—  [    cosd)  /   sin8   d6     +     sin$      /   cos8   d6    ] 
IT  o  o 


—  V     cose 

TT       S 


(III. 3. 7. 7) 


There  may  be  a  portion  of  the  cycle  when  v  exceeds  v   and  at  that 

s  K 

time  diodes  CB  and  BA  may  conduct  whilst  CD  and  DA  are  switched  off  or 
vice  versa  depending  upon  the  sign  of  v^.  The  treatment  is  also 
neglecting  non-linearities  in  the  diode  characteristic. 
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Fig. (35)   (a)  Nonlinear  characteristic  of  a  diode 
(b)  Modified  circuit  for  diode  detector 
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III. 3. 8  Diode  Detector 


The  idea  of  amplification  of  diode  detectoi"  comes  from  linear 
part  of  the  diode  characteristic  of  the  current-voltage  relation  in 
a  semiconductor. 

The  input  signal  V,  and  output  signal  V   characteristic  of  a 

diode  is  shown  in  Fig. (35a).  In  the  case  of  small  V.,  V   the 

°  1    o 

amplification  V  /V.  is  reduced.  This  causes  distortion  on  small  spin 

0    1 

echoes  and  gives  a  smaller  T„  than  the  true  value.  The  modification 
of  the  diode  detector  is  given  in  Fig.  (35b). 

For  the  system  used  in  this  dissertation  study,  an  additional 
0.3  volt  reverse  bias  is  applied  to  push  the  characteristic  toward 
the  left  so  that  a  linear  amplification  is  obtained. 

Since  the  diode  has  rectification  property,  the  information 
collected  gives  only  the  magnitude.  Thus  the  signal  averager 
(Fabritek  1062)  only  accumulates  positive  input  and  fails  to  average 
out  random  noise.  In  Fig. (36),  the  basic  structure  of  a  diode  detector 
is  shown.  The  voltage  V  is  used  to  shift  the  signal  level. 


OUTPUT 


V  DRIVES 


—►I ^~H 1      AT  POINT  kjSh-JSh 


AT  POINT  B 


SHADED  REGION  LEFT  AFTER 
DIODE  DETECTOR 


Fig. (36)   Diode  detection 
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A  diode  detector  is  also  inferior  in  tuning.  The  resonant 
frequency  is  determined  by  mixing  a  reference  signal  (from  the  pulse 
generator)  through  an  attenuator  with  the  nuclear  signal.  The  attenuator 
is  used  to  shift  the  level  of  the  signal  to  reduce  the  base  line  offset. 
When  tuned  off  resonance,  there  are  beats  on  the  free  induction  decay. 
One  can  tune  either  field  or  rf  frequency  to  remove  this  ripple  thereby 
ensuring  that  the  resonance  condition  is  fulfilled.  When  making  the 
measurements,  reference  signal  is  removed. 

In  diode  detection,  the  requirement  of  stabilities  for  field 
and  rf  source  are  not  as  critical.  Spin-spin  relaxation  times  obtained 
from  diode  detection  was  generally  smaller  by  about  30  %  than  that 
obtained  from  Gill-Meiboom-Carr-Purcell[ 36]  pulse  sequence  phase 
detection. 


III. 4   Continuous  Wave  Method 

A  Varian  DP-60  Wide  Line  spectrometer  was  used  to  record  the 
narrow  line  of  the  liquid  crystals  and  microemulsions .  The  wide  line 
of  the  liquid  crystals  contributed  by  the  deuteron  quadrupolar 
splitting  of  the  bound  D„0  was  also  recorded  by  saturating  the 
central  narrow  line. 


III. 4.1  Power  Saturation  and  Line  Width  Broadening 

Figure(37)  shows  the  relation  of  signal  to  H  (magnitude  of  rf 

field)  and  H  (magnitude  of  audio  modulation  field).  We  will  see  diagram 
m 

(a)  and  (b)  being  both  similar.  At  low  H.,  and  H  values  the  signal 
°  1      m 

height  is  proportional  to  H  and  H  ,  but  at  high  values  saturation 
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occurs  followed  by  a  decrease  in  signal  amplitude. The  corresponding  line 

width  is  no  longer  a  constant,  being  broadened  by  instrumental  reasons. 

To  avoid  such  saturation  problem,  one  should  check  the  magnitude  of 

saturation  field  before  running  each  sample  and  set  both  rf  and  audio 

field  magnitude  by  50  %  below.  When  phase  detected  at  the  modulation 

frequency,  the  familiar  derivative  spectrum  is  observed.  This  applies 

to  the  frequency  0J   is  less  than  the  linewidth  expressed  in  frequency 

unit. 

In  the  severe  saturation  a  series  of  sidebands  are  produced  at 

W   ±o),0)  ±2oj  etc.  [  Fig.  (38)]  .  The  effective  H,  at  the  sideband 
o   m   o   m  1 

frequencies  is  affected  by  H  and  can  be  expressed  as  H,  multiplied 

m  1 

by  a  scaling  factor.  The  scaling  factor  is  a  Bessel  function  of  parameter 

YHm 
(  $  =  )[37J.  Figure  (3q)  shows  that  increasing  H   requires  higher  H 

OJ  ml 

m 

to  obtain  the  same  signal  height.  In  other  words,  smaller  signal  height 

is  obtained  if  H  is  increased  and  H.,  is  kept  fixed.  (Note:  This  is 
m  1 

true  for  the  central  band  but  not  for  the  sidebands,  which  will  be 
discussed  in  next  section. ) 

For  the  first  sideband  detected  at  the  modulation  frequency 


Heffective=  ^  j2 (  p) 

rr  YH 

u     8  ■      ,1    r /    n   \        n  1   „effective  tT  .    '  m  .2 

when   M  is  small:   J(  g  )  =  g   and   H         =  H  ( ) 

-'-  -L   CO  „ 


Change  of  H   is  basically  equal  to  change  of  H  .  In  practice, 

H  is  selected  at  some  suitable  (linear  region)  value  and  final 
m 

adjustment  are  made  using  H..  . 

Thus  the  frequency  of  audio  modulation  f   determines  whether  the 
line  shape  is  in  derivative  or  absorption  mode.  For  the  convenience 
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of  comparison,  we  can  convert  the  spectral  line  width  (half  height 
width  of  spectral  line)  to  corresponding   Af .  When  phase  detected  at 
the  field  modulation  frequency,  the  familiar  derivative  spectrum  is 
obtained.  The  discussion  thus  far  applies  to  the  situation  where  the 

modulation  frequency  co  is  less  than  the  linewidth  expressed  in 

n       m 

frequency  units.  In  the  phase  angle  $   between  the  audio  output  and 
reference  is  adjusted  to  j   ,  the  largest  central  peak  is  expected 
[graph(a)  of  Fig. (40)]  .  The  change  of  polarity  of  the  central  peak  in 
this  graph  is  not  a  'derivative'  case;  this  must  be  noted  during  data 
analysing.  Graph(b)  of  Fig.  (40)  shows  the  correct  audio  phasing  at  cj>=Q 
removes  the  central  narrow  line  and  obtains  prominant  side  bands  if 
the  effective  H   is  large  enough.  Under  this  condition  the  first 

sidebands  are  the  undifferentiated  absorption  spectra,  which  are 

e  f f e  c tive 
subject  to  saturation  by  H,         in  the  usual  manner.  The  line 

width  and  signal  height  of  each  side-band  depend  on  different 

characteristic  curves  shown  in  Fig.  (41). 

Suppose  spectrum  is  more  complex  as  shown  in  graph(a)  of  Fig. (42). 

If  f   is  much  larger  either  Af,  or  Af  „  ,  we  expect  to  see  both  line 
m  1      Z 

shapes  are  in  absorption  mode  as  shown  in  graph (b).  If  f   is  smaller 

than  Af9,  we  will  see  both  in  derivative  mode  as  shown  in  graph  (c). 

If  f   is  between  Af,  and  Af„,  then  the  narrow  line  is  in  absorption 
m  1       z 

mode  and  the  wide  line  is  in  derivative  mode;  the  spectrum  is  expected 
as  shown  in  graph (d),  and  the  frequency  between  the  two  narrow7  lines 
is  still  2f  .  Graph  (e)  shows  that:  misphasing  of  the  third  case  can 
cause  other  sidebands  besides  the  first  sidebands  and  each  of  them 
with  same  width  in  the  case  of  no  saturation. 
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(a.)      £     >Af,      ?-? 
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Fig.    (40)  Line   shapes   of    the   case    :    f    >Af 
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(  a  )   TRUE  LINE  SHAPE 


(  b  )   CASE  1 

f  »Af 


2  f 


(  c  )   CASE  2 

f  <Af 


(  d  )   CASE  3 


Af2<  fm<Afl 


(  e  )   MISPHASING 
OF  CASE  3 


Fig. (42)   Line  shapes  corresponding  to  different  modulation  frequency 
f   of  signal  composing  both  narrow  and  wide  line 


85 


III. 4. 2   Rate  of  Sweep  and  Magnetic  Field  Inhomogeneity 


The  line  shap'.-  also  depends  upon  the  time  taken  to  sweep  through 

dH 
dt 


resonance.  We  may  use  AT  =  Ah/ (  j —  )  t0  characterize  this  time,  where 


dH 
AH  is  the  line  width  in  magnetic  field  units  and  ~^—   the  sweep  speed. 

at 

Figure (4 3)  shows  how  th   sweep  rate  affects  the  line  shape.  In  case (I) , 
we  obtain  normal  derivative  spectra  when  AT  is  much  larger  than  the 
time  constant,  and  distortion  of  spectra  when  AT  is  close  to  the  time 
constant.  In  case  (II),  we  obtain  further  distortion  of  spectra  when 
AT  is  smaller  than  the  time  constant.  The  symptom  will  also  be  observed 
if  AT  is  smaller  than  the  spin  lattice  relaxation  time  T-,  .  Thus  the 
sweep  rate  has  to  be  sufficiently  slow  that  AT  is  slightly  in  excess 
of  the  time  constant  or  T  whichever  is  the  larger,  whereas  the  noise 
is  inversely  proportional  to  the  square  root  of  the  time  constant. 
For  optimum  performance  the  time  constant  should  not  be  less  than  T  . 
Noise  can  be  reduced  by  increasing  the  time  constant  and  scanning  more 
slowly  or  by  making  the  two  times  comparable,  scanning  as  rapidly  as 
possible,  and  signal  averaging. 

Complication  also  arises  from  magnetic  field  inhomogeneity. 
Figure(44)  shows  that  distortion  due  to  magnetic  field  inhomogeneity 
may  have  a  similar  symptom  e.g..  if  the  field  has  a  quadrupole  component 
which  produces  an  asymmetric  line  shape.  We  can  distinguish  this  by 
checking  the  scan  in  both  directions.  If  the  polarity  of  the  spectra 
does  not  change, this  is  due  to  magnetic  field  inhomogeneity.  If  only 
polarity  change  and  distortion  is  in  the  same  side  as  before,  this  is 
due  to  rapid  passage. 

The  CW  measurements  recorded  in  this  dissertation  employed  only 
slow  passage  conditions.  Though  the  fast  passage  condition  is  not 
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CASE  (I)   AT  >  TIME  CONSTANT 


(a)   AT  »  TIME 

CONSTANT 


(b)   AT  X  TIME 

CONSTANT 


CASE  (II)  AT < TIME  CONSTANT 


(  SAME  SYMPTOM,  IF  AT<  T   WHERE  T   IS  SPIN-LATTICE 

RELAXATION  TIME) 


Fig. (43)   The  effect  of  slow  and  fast  sweep  upon  the  spectri 
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CASE  A   DUE  TO  FIELD  INHOMOGENEITY 


—  H  -> 


<-H 


TIME   > 


TIME 


CASE   B       DUE  TO  RAPID   PASSAGE 


—  H     ~> 


TIME 


TIME 


Fig. (44)   The  distinction  of  line  distortion  due  to  magnetic  field 
inhomogeneity  and  rapid  passage 


strictly  relevant,  we  still  discuss  it  a  little.  This  is  the  case 
when  AT  is  smaller  than  T  .  Under  these  conditions  the  initial  signal 
height  is  proportional  to  the  amplitude  of  the  rf  field.  Before  a  new 
equilibrium  (i.e.  saturation)  can  be  established  the  field  has  passed 
through  the  resonance  condition.  In  practice  this  condition  cannot  be 
completely  achieved,  the  later  parts  of  the  signal  will  be  saturated 
and  appear  with  reduced  value.  In  this  circumstance,  fast  passage  will 
increase  the  linear  region  of  signal  height  versus  H   characteristic 
and  allow  the  observation  of  a  bigger  signal  than  would  otherwise  be 
possible  but  it  will  have  distorted  line  a  shape  as  shown  in  case  (II) 
of  Fig. (44). 


Ill . 5  Noise 

A  major  problem  encountered  by  the  NMR  spectroscopist  is  the 
extraction  of  the  signal  froir  a  noisy  spectrum.  Noise  originates  at 
each  stage,  but  obviously  that  produced  in  the  early  stages  of  a 
high  gain  amplifier  is  the  most  significant.  In  NMR  transmitter  noise 
tends  not  to  be  significant  as  (i)  in  pulse  NMR  the  signal  is  observed 
only  when  the  transmitter  is  disconnected  (ii)  in  CW-NMR  the  bridge, 
mechanical  or  electrical,  attenuates  the  transmitter  power  and  noise 
coupled  into  the  receiver.  There  are  usually  some  tuning  elements  in 

the  receiver .which  give  it  a  bandwidth  Af  .  The  noise  is  proportional 

1/2 
to  (Af,)    ,  and  in  principle  could  be  reduced  by  reducing  Af  ,  but  as 

this  would  raise  problems  in  maintaining  the  alignment  of  the  receiver, 

Af   is  usually  made  large.  The  output  of  the  radio  receiver  is  then 

detected  by  either  diode  or  phase  detector.  The  effect  is  to  mix  the 
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output  of  the  radio  receiver,  signal  noise,  with  a  reference  signal 
at  frequency  f   or  with  a  leakage  signal  at  frequency  f  .  Application 
to  a  non-linear  stage  results  in  the  generation  of  sum  and  difference 
frequencies  and  multiples  thereof.  The  components  near  zero  frequency 
are  detected.  This  procedure  results  in  a  signal  at  zero  frequency, 
in  the  absence  of  field  modulation  or  at  the  field  modulation 
frequency  in  it's  presence,  together  with  noise  up  to  a  frequency  of 
Af1 .  This  assumes  that  the  additional  noise  introduced  in  the  detector 
stage  is  negligible.  In  pulsed  NMR  this  signal  is  then  displayed,  a 
low  pass  filter  with  cut-off  at  Af~  can  be  inserted  to  limit  the 
bandwidth  and  thus  the  noise.  To  prevent  signal  distortion  Af   should 
exceed  the  highest  significant  frequency  components  of  the  Fourier 

Transform  of  the  signal,  i.e.  if  the  characteristic  time  of  the  free 

*         1 

induction  decay  is  T„  ,  Af„>  —■ *  . 

In  CW-NMR  field  modulation  is  usually  employed  and  in  order  to 

eliminate  changes  in  receiver  output  arising  from  changes  in  transmitter 

power  bridge  balance,  amplifier  gain  etc.  use  is  made  of  the  fact 

that  the  signal  depends  upon  magnetic  field  whereas  these  other  factors 

do  not,  and  that  a  sinusoidal  change  in  magnetic  field  at  frequency  f 

m 

will  produce  a  sinusoidal  change  in  NMR  signal  at  frequency  f  .  Mixing 

with  a  reference  signal  at  frequency  f   again  generates  sum  and 

difference  frequencies  and  the  zero  and  near  zero  frequency  components 

can  be  selected  using  a  suitable  phase  detector.  It  is  then  convenient 

to  limit  the  bandwidth  Af   in  the  subsequent  d.c.  stage.  Again  care 

must  be  taken  to  ensure  that  Af_  is  sufficiently  large  not  to  distort 


the  lineshape.  If  the  sharpest  spectral  feature  has  a  width  AH  and  the 

3' 


sweep  rate  is  dH/dt, then  Af  >  (AH)   (dH/dt) 
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Because  of  problems  of  receiver  alignment  it  is  invariably  more 
convenient  to  limit  the  bandwidth  at  the  receiver  stage.  Sources  of 
noise  are  (i)  Johnson, (ii)  Shot, (iii)  Partition, (iv)  Flicker.  The 
first  three  are  white  broad  band  noise,  whereas  the  last  comprises 
low  frequency  components.  The  modulation  frequency  is  made  as  high  as 
possible  to  minimise  it's  contribution. 


TT 

III.  6  T  Measurement  with  it  "u   t  ^ ■  —  Pulse  Sequence 


III. 6.1  Disadvantage  of  Cross-Over  Method 


The  easiest  way  to  measure  T   is  by  using  the  cross-over  method, 
which  finds  the  delay  time  giving  a  null  signal  after  the  tt/2   pulse. 
If  the  magnetization  M   is  flipped  by  a  precise  tt  pulse,  we  can  adjust 
the  delay  time  t  to  give  a  null  free  induction  decay  after  the 
following  tt/2   pulse.  It  is  easy  to  show  that  T   is  equal  to  t/ln2. 
However,  in  practice  it  is  difficult  to  tune  a  precise  tilting  angle 
of  TT  and  this  affects  the  measurement  of  T-  as  shown  by  the  analysis 
given  below. 

From  the  Bloch  equation,  the  rate  of  change  of  the  z  component 
of  magnetization  is 


dM    M  -  M  dM 

z    o z  z dt  ,TTT  ,  ..  ..  s 

IT  =  -Y—  or  m^¥~=t"        (iii. 6. i.i) 

1  o    z      1 

-In [11  -  M  (t)]  = 1-  C   where  C  is  a  constant  depending  upon 


Tl 


initial  condition. 
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Mq  -  Mz(t)  =  exp(-  ~)    exp(-C) 


(III. 6. 1-2) 


Initially:  after  a  pulse  of  6  [Fig.  (45)] 


^ 


M  cos6 
to 


Fig. (45)   Magnetization  flipped  down  by  an  angle  0  from  initial  direction 


M  -  M  cos0  =  exp(-C) 
o    o 


M  -  M  (t)  =  M  (  1  -  cosG  )  exp(-  -— ) 
o     z       o  T 


(III. 6. 1.3) 


Therefore 


M  (t)  =  M  (  1  -  (  1  -  cos6  )  exp(-  — ) 


(III. 6. 1.4) 


The  cross-over  method  measures  the  time  when  M  (t)  =  0 


i.e.  0  =  M  [  1  -  (  1  -  cosG  )  exp(-  — ) 

T1 


This  gives   T..  =   :  j~. ■ jr~r 

1    ln(  1  -  cos0  ) 


(III. 6. 1.5) 


Thus  the  answer  depends  on 


If  6  =  IT   then  T  =  t/ln2 
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Suppose  the  tilt  angle  0  is  not  7T,  then  the  existing  error  is 


t.  =  t  (  ~-  -  T~rr sr)  (in. 6. 1.6) 

1       ln2    ln(  1  -  coso)  J 


Figure (46)  shows  that  two  curves  give  different  cross-over 
position,  even  though  subject  to  same  T,  relaxation  scheme. 


III. 6. 2   Ringing 

One  of  the  errors  of  measuring  the  free  induction  directly 
comes  from  ringing  of  the  pulse.  Figure (47)  shows  how  the  tail  of  the 
pulse  can  mix  with  the  free  induction  decay.  The  procedure  used  was 
to  find  the  intercept  of  the  extension  of  the  FID  with  an  "ideal" 
pulse  and  then  measure  the  amplitude  between  this  intercept  point 
and  the  base  line.  The  width  of  the  FID  of  the  same  sample  varies  at 
different  temperatures.  A  narrow  FID  with  steep  slope  enhances  the 
error.  Failure  to  compensate  will  result  in  error,  e.g., a  20  %  smaller 
T1  in  severe  case  [Fig. (48)]. 

Ringing  is  a  basic  damping  problem  of  the  LRC  circuit.  For 
simplification  of  the  problem,  one  can  denote  the  total  resistance, 
capacitance  and  inductance  as  R,  L  and  C  respectively.  The  damping 
for  t  <  0  can  be  ignored  because  of  the  large  driving  power  of  the  rf 

source.  For  t  >  0,  there  are  three  types  of  damping,  namely: 

2   2        -1 
(i)  R  /4L   >  (LC)     This  is  the  condition  for  OVER  DAMPING. 

The  voltage  has  two  mixed  exponential  decay 

schemes  and  decreases  to  zero  at  t  =  co. 

(ii)  R  /4l/  =  (LC)"1   This  is  the  condition  for  CRITICAL  DAMPING, 
the  system  is  said  to  be  DEAD  BEAT.  The 
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PULSE  (  ABOUT  20/lsec) 


INCORRECT  POINT  TO  TAKE  DATA 


TAIL  DUE  TO  PULSE  RINGING 


CORRECT  POINT  TO  TAKE  DATA 


POINT  OF  DISCONTINUITY 


NATURAL  EXTENSION  OF  FID 


TIME 


Fig. (47)      The    observed    FID   as   displayed   on   the  X-Y   recorder 
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voltage  has  only  a  single  exponential  uecay 
scheme  and  decreases  to  zero  in  the  shortest 
time.  This  is  what  we  prefer. 

(iii)  R2/4L2  <  (LC)_1  This  is  the  condition  for  DAMPED  OSCILLATION. 
and  the  case  we  are  most  concerned  with.  It 
the  R/2L  is  small,  the  voltage  takes  a  long 
time  to  damp  out,  thus  causes  serious 
ringing. 

Since  the  phase  at  which  the  radio  wave  (contained  within  the  ^./2 

pulse  envelope)  is  cut  by  the  pulsing  gate  is  beyond  control,  it  is 

difficult  to  provide  an  initial  condition  to  remove  ringing.  Therefore 

we  must  increase  the  R/2L  ratio  to  have  fast  damping.  The  resistance  R 

comes  predominantly  from  the  probe  in  most  cases,  but  other  resistance 

from  the  transmitter  and  receiver  contributes  also.  To  construct  a 

coil  with  finer  coil  wire  can  increase  resistance  dramatically,  but 

decreases  the  quality  factor  oj  L/R.  Fortunately,  the  deuteron  signal 

o 

is  strong,  and  the  signal-to-noise  ratio  in  most  of  the  temperature 
range  is  good;  hence  a  large  quality  factor  is  comparatively  not  so 
critical. 

III. 6. 3  Dead  Time 

If  a  tt/2  pulse  of  time  duration  At  is  given  to  a  resonant  sample. 
the  free  induction  decay  starts  from  At/2  in  principle.  In  practice 
the  FID  is  observed  starting  f rom(At/2)4b.  No  information  is  observable 
within  the  period  b.  This  time  interval  is  called  DEAD  TIME.  The  dead 
time  is  caused  by  the  recovery  of  the  system.  Bad  tuning  increases  the 
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dead  time  tremendously;  therefore,  one  must  accurately  tune  the  system 
to  reduce  dead  time  and  preserve  more  information.  A  well  tuned  pulse 
must  have  the  shortest  decay  and  rise  time  by  decreasing  time  constant. 
Unfortunately,  there  is  a  limit  for  doing  this.  Large  capacitance 
needs  large  inductance  of  the.  sample  coil  for  balance  to  avoid  serious 
ringing,  as  has  been  discussed  in  the  last  subsection.  To  obtain  a 
strong  signal,  we  made  the  sample  as  large  as  the  physical  dimensions 
of  cryostat  permitted;  hence  the  geometrical  size  of  the  sample  coil 
was  predetermined.  A  large  inductance  then  requires  finer  wire  if  the 
geometry  of  the  coil  is  fixed  and  this  increases  resistance  which  in 
turn  reduces  gain.  Therefore  we  can  only  tune  the  system  by  a  compromise 
among  L,  R,  and  C. 

The  signal  developed  in  the  receiver  coil  is  proportional  to  the 
Q  factor  of  the  coil  and  should  therefore  be  high.  The  dead  time  is 
inversely  proportional  to  the  bandwidth  and  therefore  the  bandwidth 
should  be  high  and  the  Q  factor  low.  Design  of  a  suitable  coil  system 
requires  a  compromise  between  these  competing  factors. 

For  the  20  ysec  pulse  of  this  experiment,  the  dead  time  is  about 
150  psec.  This  is  short  compared  with  the  relaxation  time  of  ten  to 
a  few  hundred  msec  over  most  of  the  temperature  range.  Below  -40  C, 
most  of  the  samples  studied  have  their  relaxation  time  below  2  msec, 
and  a  larger  error  is  expected. 
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7J 

III.  7   T„  Measurement  with  ---  ^  x  ^  TT   Pulse  Sequence 


IT 
The  -T-  ^  x  ^  IT  pulse  sequence  was  employed  to  measure  the  envelope 

of  the  spin  echoes.  The  formation  of  spin  echoes  has  been  mentioned 

in  the  basic  theory  in  chapter  II.  To  obtain  a  good  echo,  correct 

tuning  in  both  pulse  and  signal  is  essential.  In  this  section,  we  shall 

discuss  the  analysis  of  the  data  for  the  liquid  crystal  which  exhibits 

two  different  relaxation  times. 

Since  a  diode  detector  was  used  to  measure  the  spin-spin  relaxation 
time  1    ,    all  echoes  are  subject  to  a  basic  line  shift  (chapter  III, 
section  III. 3. 8);  the  shift  can  be  removed  by  self  adjustment  in 
computer  programming.  In  graph(a)  of  Fig. (49),  we  can  see  two  types 
of  distortion  causing  a  non-exponential  decay  curve  of  the  spin  echo 
in  the  liquid  crystals  phase.  The  distortion  in  the  tail  is  especially 
subject  to  the  base  line  shift  and  non-linearity  of  the  diode  for 
small  input  signal.  The  second  factor  is  much  improved  by  using  a  DC 
voltage  supply  of  approximate  0.3  volt  to  shift  the  input  signal 
(section  III. 3. 8).  The  distortion  near  the  origin  arises  is  because 
two  pulses  are  too  close. 

In  graph(b)  of  Fig. (49),  liquid  crystal,  a  typical  example,  is 
shown  for  the  hexanol  W/0  =1.2  with  S/0  =  0.2.  There  are  two  features, 
one  for  tightly  bound  water  of  the  same  kind  as  shown  in  graph  (a), 
the  other  due  to  isotropic  water.  Therefore  data  taken  for  least 
square  fit  must  be  treated  separately.  When  we  take  height  of  echoes, 
the  distortion  area  of  graph (b)  must  be  excluded,  otherwise  larger  T 
for  tightly  bound  water  and  smaller  T„  for  isotropic  water  are  expected. 


DISTORTION  DUE  TO  VERY 
CLOSED  PULSE  SEQUENCE 
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(a)  HEXANOL  LIQUID  CRYSTAL 
W/0  =  1.4  ;  S/0  =  0.2 


ERROR  DUE  TO  BASE  LINE  SHIFT  AND 
NONLINEARITY  OF  DIODE  DETECTOR 


0.5  ms 


(b)    HEXANOL  LIQUID   CRYSTAL 
W/0  =   1.2    ;    S/0  =   0.2 
t  =   12.1°C 


0.5  ms 


Fig. (49)    Envelope  of  spin  echoes  of  two  typical  liquid  crystals 
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III. 8  Least  Square  Fit 

The  amplitude  of  the  free  induction  decay  following  a 

pulse  is  given  by  y  =  A  -  B  exp (-  —  )  .  T  was  deduced  from  a  least 

1 
square  fit  of  this  equation  to  experimental  data  following  a  procedure 

described  by  Sass  and  Ziessow[38]. 

The  least  square  fit  is  based  on  the  assumption  that  the  error 

of  the  measured  signal  has  a  Gaussian  distribution.  The  error  square 

sum[39][40]  is  given  by 

N 
E  =L     [  y.  -  A  -  B  exp(  at.)]  =  0  (III. 8.1) 

i=l   X  X 

where  v.  is  the  ith  measurement  of  the  magnetization  at  delay  time  t., 

J  ±  __  o  j  x 

and   a  =   -  —  .    To  minimize.  E,    the   three  partial   derivatives  ttt1   j      T,T   > 
T  '  dA  dB 

,    8E 
and  t—     are   set   to   zero. 
9a 


1     ^F  N 

-fff  =  £     [   y±   -  A  -  B  exp(  at.)]    =   0  (III.  8. 2) 

i=l 

-t!t  =  ^     [    y.    -  A  -  B   exp(   at.)]    exp(    at.)    =    0  (III. 8. 3) 

L      OB  .  _,  X  1  1 

1       3K  N 

-  f^TT  -    Z     [   Y.    -  A  -  B   exp(  at.)]    t.    exp  (  at.)   =   0    (III.  8. 4) 

Lb    da  .     ,  1  11  1 

1=1 


One  then  solves  for  equations  (III. 8. 2)  and  (III. 8. 3)  as  a 
function  of   a.  Introducing  the  abbreviations 


N 

2,  x.  =   {x}      e.  =  exp (  at.) 

i=l 

t.*    ■              a         (y)         {y){e}    -   N{ye}    ,    ,  . 

one  obtains        A  =  — *: *-= = l „{e}  (III. 8. 5) 

N  UZ{eZ}    -  N{e}2 
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on 


{y}{e>  -  Nfre}  (m<8>6) 

N{e^}   -    {e}2 


Substitute  A  and   B    in   equation    (III. 8. 4)    yields 


G  =    (  N{e2}   -    {e}2){yte}  +    ({ye}{e}    -    {y}{e2}){te} 

+    ({y}(e}   -  N{ye}){te2}   =    0  (III. 8. 7) 


where  G   is   a  refined   value   of      as 

The  null   position  of   G   is    iteratively   determined  with    the   Newt 
method : 

a(k+l)   «   aUO    -    [    G(k)/G'(k)],  (III. 8. 8) 

where  the  index  k  indicates  the  k  th  choice  of  a  and  G'  is  the  partial 

.   .    .    3G 
derivative  tt-  ! 

da 

G'  =  (  N{te2}  -  {te}{e}){yte}  +  (  N{e2}{yt2e} 

+  ({yeHte}  -  {y}{te2}){te}  +  ({ye}{e}  -  {y}{e2 » {t2e} 
+  2  ({y}{e}  -  N{ye}){t2e2}  .  (III. 8. 9) 

The  iteration  (III. 8. 8)  was  programmed  in  BASIC  for  a  PDP-8/A 
CLASSIC  computer.  Input  data  are   the  number  of  signals  N,  the  free 
induction  decay  amplitude  y.,  the  delay  time  t.  and  a  guessed  initial 
value  for  T  .  If  the  choice  of  T   is  within  ±50  %  of  the  final  T   value, 
no  more  than  5  iteration  steps  are  usually  required  to  fit  T, ,  A  and 
B  with  an  accuracy  of  ±1  %. 

Two  programs  are  attached  in  the  appendix.  The  second  program  is 
the  same  as  the  first  one  except  that  provision  is  made  for  entering 
data  from  a  large  number  of  measurements  at  one  time  rather  than  one 
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set  separately  for  each  fit.  The  program  self adjusts  the  base  line 
and  does  not  need  to  know  y  at  t  =  °°.  For  fast  convergence  the  data 
should  include  points  in  the  region  of  fast  change  of  signal,  i.e. 
roughly  from  t  =  0  to  I,  .  The  correspondence  of  notation  between  the 
analysis  above  and  the  program  is  given  by  the  following  equivalences 


y.  E  A(I),     t.  E  T(I),     {e}  E  S  (0) ,     {e2}  E  S(l) 


{y}  E  S(2),     {ye}  E  S(3),     {yte}  E  S(4),     {te}  E  S(5), 

{te2}  E  S(6),     {yet2}  E  S(7),     {et2}  E  S(8),     {e2t2}  =  S(9), 

A  E  H,    BEM,    a  =  G,    G  =  Gl,    G '  E  G2 

The  program  assumes  an  exact  determination  of  t . .  The  error  in  t .  is 
r   o  x  ^ 

within  1  x  10   sec;  therefore,  the  influence  on  y.  can  be  negligible. 
The  least  square  fit  method  is  independent  of  tilt  angle  6. 

The  Wang  pulse  programmer  used  in  this  study  gives  an  unstable 
delay  time  for  t.  larger  than  216.15  ms .  This  causes  slow  converge  and 
errors  in  the  case  T1  is  larger  than  200  ms.  Therefore  the  semilogarithmic 
graphical  method  was  employed  also  by  plotting  [  y  (°°)  -  y.(t)]  on  a 
semilogarithmic  scale  against  delay  time  and  calculating  T   from  the 
slope.  This  method  is  also  free  from  the  inaccurate  adjustment  of 
the  tilted  angle  G,  but  needs  to  know  y.  at  t  =  »,  which  equals  y   at 
t  =  0. 

The  same  computer  program  was  used  for  the  least  square  fit  of 
the  spin  envelope  to  calculate  T„  values.  The  program  selfadjusts  the 
base  line  and  does  not  need  to  know  the  very  small  echo  amplitude  for 
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the  very  large  delay  time.  The  semilogarithmic  method  was  also  employed 
for  T„  larger  than  2  50  ms  for  the  same,  reason  mentioned  in  the  last 
paragraph. 


III.  9  Visual  Observation 

The  purpose  of  the  visual  observations  was  to  check  birefringent 
phenomena  and  phase  separation  in  order  to  compare  with  the  relaxation 
time  measurements,  which  were  performed  in  a  metal  cryostat  where  the 
samples  were  not  visible.  A  limitation  of  the  visual  observations 
occurs  at  lower  temperature.  Below  the  melting  point  of  hexadecane  18  C. 
the  samples  become  an  opaque  milk-like  liquid  and  visual  observation 
is  no  longer  possible. 

The  apparatus  used  is  very  simple,  as  is  shown  in  Fig. (50).  The 
system  is  basically  a  gas  flow  cryostat  in  which  gas  at  a  regulated 
temperature  is  passed  through  a  glass  U  tube  containing  the  sample 
inside  a  copper  can.  For  observation  above  room  temperature,  we  passed 
air  through  a  high  power  electrical  heater  and  then  through  the. 
cryostat.  Both  flow  rate  and  current  of  the  heater  were  controlled  to 
raise  the  temperature  of  the  copper  cryostat.  For  the  observation 
below  room  temperature,  we  introduced  cold  nitrogen  gas  flow  by 
controlling  the  current  is  a  resistor  submerged  in  an  external  liquid 
nitrogen  bath.  The  pyrex  glass  U  tube  was  mounted  inside  a  large  glass 
dewer.  Diagram  of  the  copper  cryostat  is  shown  in  Fig. (51).  It  has  an 
external  diameter  of  one  inch  and  internal  one  of  0.5  inch  to  contain 
the  sample.  Glass  slits  with  a  width  1/8  inch  located  diametrically 
opposite  on  the  cylinder  provide  a  light  path  for  viewing  the  sample. 
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Two  polarized  plates  were  located  in  opposite  sides  of  the  slits. 

The  copper  cap  limits  heat  radiation.  Two  pairs  of  thermocouples 
are  located  inside  to  measure  temperature  of  the  sample  at  different 
positions.  This  cryostat  gives  a  homogeneous  temperature  distribution 
all  over  the  sample  because  of  the  excellent  thermal  conductivity  of 
copper.  The  axial  temperature  gradient  can  be  reduced  if  the  time  for 
heating  or  cooling  is  sufficiently  long.  In  general,  it  takes  5-10 
minutes  to  reach  thermal  uniformity  near  room  temperature  and  25-40 
minutes  at  70  degrees  from  room  temperature.  The  glass  dewer  in  Fig. (50) 
is  sealed  with  tape  and  a  vapor  absorber  material  is  placed  on  the 
bottom  to  prevent  ice  crystalization  on  the  surface  of  the  glass  U 
tube  when  the  temperature  is  reduced.  The  thermocouples  were  connected 
to  a  potentiometer  by  a  DPDT  switch.  An  ice  water  mixture  was  used  as 
reference  junction. 

In  the  case  of  the  microemulsions ,  light  is  scattered  isotropically 
by  the  small  water  droplets;  hence  we  expect  to  see  uniform  illumination 
as  the  polarized  plates  are  rotated  and  to  have  maximum  intensity  when 
the  plates  are  orthogonal  to  each  other. 

In  the  case  of  the  liquid  crystals,  the  light  is  reflected 
anisotropically  because  of  their  grating-like  structure.  When  the 
path  difference  of  light  is  an  intefral  multiple  of  the  wavelength, 
there  is  interference  and  we  observe  colorful  birefringence.  Since 
the  gap  between  two  layers  is  of  the  order  10  A,  which  is  small 
compared  with  the  wavelength  of  visible  light,  the  expected  angular 
distribution  is  close  to  90  .  If  we  observe  a  completely  dark  sample 
with  two  plates  orthogonal  to  each  other,  this  means  that  the  gap 
between  layers  is  so  narrow  that  the  light  beam  is  not  transmitted. 
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This  special  case  has  been  observed  for  a  pentanol  sample  of  W/0=0.6 
with  S/0=0.2. 

After  phase  separation  at  high  temperature,  the  liquid  crystal 
takes  2  to  10  hours  to  return  to  its  original  structure.  If  the  sample 
is  taken  out  of  the  cryostat  and  shaken  well,  the  structure  is  restored 
immediately.  This  is  consistant  with  the  experience  that  T  measurement 
of  the  liquid  crystals  in  the  phase  separation  region  at  high 
temperature  shows  thermo-hysterisis . 

The  visual  observations  were  carried  out  from  room  temperature 
down  to  the  freezing  point  of  hexadecane(18  C)  and  from  room  temperature 
up  to  100°C. 
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Fig. (51)     Cryostat  for  visual  observation 


CHAPTER  IV  RESULTS  AND  CONCLUSIONS 

I V . 1  The  Micro emulsions 

IV. 1 . 1  The  Validity  of  Bound  Water  Model 

The  bound  water  model  proposed  has  been  discussed  in  section  II . 4 
of  chapter  II.  Figure (52)  shows  the  relationship  between  line  width 
and  radius  of  the  microemulsion  spheres.  The  volume  ratio  of  oil/water 
(0/W)  can  be  converted  as  in  the  following: 


I-  (f  )(|)  =  (|rl(|)  (iv.i.i.D 

3 
where  S/0  is  the  ratio  of  surfactant  to  oil  in  gm/cm  ,  and  has  the 

value  0.2. 

If  there  are  n  water  droplets,  the  total  water  volume  will  be 

4irr 
n  — - —  ,  where,  r  is  the  radius  of  each  spherical  droplet.  Since 

potassium  oleate  has  a  highly  polar  terminal  group  and  a  hydrophobic 

alkyl  chain,  it  is  localized  on  the  interfacial  surface  of  the  water 

2 
droplet.  Its  weight  will  be  nM  4irr  C/A  ,  where  M  is  the  molecular 

weight,  C  the  fraction  of  the  total  area  occupied  by   the  potassium 

oleate,  and  A   the  effective  area  occupied  by  each  potassium  oleate. 

The  reason  for  introducing  this  is  because,  the  localized  cosurfactant 

hexanol  may  occupy  the  remainder  of  the  area. 

Therefore,   S    (  ^^2^  /A  )/(n^r3/3)  =  3M  C/A  r     (IV. 1.1. 2) 
w  s   s  s   s 
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So,  the  ratio 

15  M  C 

H.  =  -2 —  .  i  (IV.  1.1.  3) 

W     A     r 
s 

-1 

is  proportional  to  r 

Since  simple  exponential  relaxation  is  observed  for  water  in  all 
the  microemulsion  samples,  exchange  between  this  "bound,   or  motionally 
restricted  layer  at  the  interfacial  film  and  the  "free"  (or  bulk  liquid- 
like) water  in  the  droplet  interior  must  be  rapid[41].  In  this  case, 
the  observed  water  relaxation  time  will  be  a  weighted  average  of  the 
relaxation  time  of  "free"  and  "bound"  water  as  given  by 

<  ~>  ,         =  P.    ~  +    (  1   -  P,     )   ~  (IV. 1.1. 4) 

T2      obs  b   Tb  T 

2 

Substituting     P,     -  -/r      ,-  (IV.  1.1.5) 

b  4  3 

3TT  r 

where  Ar  is  the  thickness  of  bound  water,  one  obtains 

<  ~-  >  .        -  K+   3  Ar(  —  -  ~  )  -  (IV. 1.1. 6) 
T   obs.     f          b    f   r 

2         2    2 

A  linear  relationship  between  observed  relaxation  rate  (or  line- 
width)  and  the  reciprocal  of  the  drop  radius  is  therefore  predicted. 
This  equation  also  predicts  that  if  a  range  of  drop  radii  occurred  in  a 
given  sample  at  a  fixed  temperature  there  would  be  a  range  of  T„  values. 
The  relaxation  data,  however,  provide  no  evidence  of  multicomponent 
relaxation  and  thus  indicate  that  the  drops  are  of  rather  uniform  size 
for  fixed  composition  and  temperature.  This  has  been  confirmed  by  a  light 
scattering  experiment [8] .  If  the  exchange  of  water  molecules  between 
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droplets  is  slow,  non-exponential  relaxation  would  be  predicted.  The 

slopes  of  Fig. (52)  depend  upon  both  — —  and  Ar  .  If  either  is  known 

T 

2 
independently,  the  other  could  be  calculated. 

Mathematically,  the  relationship (IV. 1. 1. 5 )  is  no  longer  true  when 

r  is  comparable  with  Ar .  If  the  outer  radius  of  the  droplet  is  r   and 

the  radius  of  the  bulk  water  core  is  r  ,  then  the  probability  to  be 

"bound"  in  the  limiting  case  will  be 

4_  ,  3    3* 

lim  P  =  3  UU1   V   =  1    (  r  -»■  Ar  ) 

ro+0        4   3  l 

2         3  ^ 

Since  the  quantity  of  the  bulk  water  then  becomes  negligible, 

the  second  term  of  f ormula(IV. 1. 1. 4)  vanishes.  Therefore  <  =—  >  ,   =  — r- 

T   obs.    b 

2 
When  we  plot.  J.og  <T„>  ,    against  radius  of  the  droplet  at  fixed 
2   obs. 

temperature  for  different  microemulsion  samples,  we  can  extend  the  line 
(fairly  linear)  to  intercept  the  ordinate  to  obtain  T  . 

From  Fig. (52),  we  can  see  that  graphs  of  the  line  width  against 
0/W  ratios  at  fixed  temperature  are  linear.  If  we  extend  this  line, 
the  interception  at  r  =0,  should  be  the  same  as  the  intrinsic  line 
width  of  DO.  Those  data  obtained  by  the  pulsed  method  are  calculated 
as  the  full  width  at  half  height  =  (tt  T  )   ,  by  assuming  a  Lorentzian 
shape.  As  the  observed  deuteron  spin-spin  relaxation  was  single  component 
exponential  this  seems  to  be  a  valid  assumption.  Linewidth  data  from 
the  continuous  wave  measurement  are  somewhat  larger,  but  exhibit  the 
same  linear  relation.  The  most  likely  explanation  for  this  difference 
in  linewidth  between  the  CW  and  pulse  methods  is  that  the  CW  measurements 
are  broadened  by  magnetic  field  inhomogeneity ,  whereas  the  Carr-Purcell 
pulse  sequence  compensates  for  field  inhomogeneity.  Another  factor 
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influencing  the  accuracy  of  the  CW  measurements  was  the  temperature 

control.  The  gas  flow  type  controlling  system  only  gave  a  fairly 

uniform  temperature  in  a  small  region  of  about  1  cm,  and  therefore  a 

thermal  gradient  exists  inside  the  sample. 

The  radii  of  different  microemulsions  are  listed  in  Table(2)  in 

chapter  I.  From  the  known  T   ,  T   (obtained  from  intrinsic  bulk  DO 

sample),  <  T„  >  .    ,  and  r,  we  may  obtain  Ar-  1  A  for  each  sample. 
2   obs. 

The  calculated  result  is  listed  in  Table(5).  T   is  8,  12,  20,  and  50  ins 
at  5,  15,  28,  and  57°C,  which  are  plotted  in  Fig. (53) .  Thus,  the 
water  relaxation  time  data  suggest  that  on  the  average,  about  1  A 
thick  layer  of  water  molecules  in  the.  water  droplet  is  restricted  in 
motion  by  interaction  with  the  surfactant  polar  groups.  This  corresponds 
roughly  to  a  "monolayer"  of  water.  This  is  also  true  over  a  wide  range 
of  temperature.  The  apparent  deviation  of  the  sample  W/0  =  0.1  can 
be  explained  that  the  composition  was  affected  by  the  freeze-pump-thaw 
technique  in  such  a  case  of  small  water  content. 

The  experimental  data  strongly  support  the  simple  geometrical 
bound  water  model  for  microemulsions. 

Table(5)   The  thickness (in  A)  of  the  motional  restricted  layer  of  the 
microemulsions  at  different  .  temperatures  and  various  radii 

q  n  n  n  r* 

r  in  A  at 

W/ 0=0.1  (r  =  16  X  ) 

0.2  (r  -  23  X  ) 

0.4*(r  =  28  X  ) 

0.6  (r  =  36  X  ) 
*  A  reasonable  guess  from  other  sizes  calculated  by  L. Prince  and  J.Schulman 
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1.80 

2.03 

1.98 
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1.06 

0.98 

1.10 

1.15 
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0.86 

1.02 

1.10 

0.90 
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1.04 
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Fig.  (52)   The  relation  between  line  width  and  radius  of  the  microemulsions 
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IV . 1 . 2   The  Rotational  Diffusion  and  Chemical  Exchange 

The  T  measurement  as  a  function  of  temperature  for  pure  D„0(99.8%) 
compared  with  the  typical  microemulsion  W/0  =   0.1  is  shown  in  Fig.  (54). 
The  calibration  data  (indicated  by  ■ )  is  about  the  same  as  that  obtained 
by  Hindman[25].  Below  the  freezing  point,  the  pure  DO  sample  does 
not  provide  a  visible  signal,  because  quadrupolar  and  dipolar  interactions 
in  the  rigid  lattice  of  D„0  give  a  broad  line  on  the  order  of  250  kHz 
wide.  The  microemulsion  curve  is  generally  parallel  to  the  calibration 
curve  except  in  the  wide  range  of  supercooling  from  -4  C  to  -44  C. 
Supercooling  is  expected  to  be  prevalent  because  in  bulk  water  freezing 
is  normally  induced  by  heteronuclei.  In  the  microemulsions  the  water 
is  distributed  in  the  form  of  droplets,  most  of  which  will  not  possess 
a  suitable  impurity.  The  parallelness  of  the  relaxation  behavior  implies 
that  the  core  water  in  the  droplets  also  experience  the  same  two 
processes  as  demonstrated  in  bulk  water[25],  namely,  (I)  hydrogen  bond 
breaking  process  at  lower  temperature  and  (II)  Brownian  rotational 
diffusion  process.  The  difference  in  relaxation  occurs  because,  motion 
of  molecules  in  the  bound  layer  is  restricted  even  though  the  exchange 
between  bound  and  free  molecules  is  rapid.  We  assume  that  the  influence 
of  paramagnetic  impurities  is  negligible.  Let  us  look  at  the  T  behavior 
of  four  microemulsions  with  W/0  volume  ratio  0.1,  0.2,  0.4,  0.6  shown 
in  Fig. (55),  Fig. (56),  Fig. (57)  and  Fig.  (58).  All  values  of  T2  above 
40  C  have  tendencies  to  a  decreased  temperature  gradient;  this  is 
caused  by  rotational  motion.  When  the  temperature  increases,  the  thermal 
expansion  increases  the  intermolecular  distance",  the  rotational  barrier 
(indicated  by  activation  energy)  becomes  smaller [Fig . (16)  ] ;  therefore, 
the  rapid  rotational  motion  dominates.  The  thermal  rotational  diffusion 
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Fig. (55)   The  temperature  dependence  of  relaxation  times  of  the 
hexanol  microemulsion  W/O  =0.1  with  S/O  =0.2 
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Fig. (56)   The  temperature  dependence  of  relaxation  times  of  the 
hexanol  microemulsion  W/0  =0.2  with  S/0  =0.2 
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Fig. (57)   The  temperature  dependence  of  relaxation  times  of  the 
hexanol  microemulsion  W/O  =0.4  with  S/O  =0.2 
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Fig. (58)   The  temperature  dependence  of  relaxation  times  of  the 
hexanol  microemulsion  W/0  =0.6  with  S/0  =0.2 
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is  a  continuous  T  decreasing  process  apparently  as  compared  with 
f igure(16) . 

An  anomaly  occurs  in  the  T  data  of  hcxanol  microcmulsions  in  the 
region  of  50  C.  The  most  probable  explanation  for  this  rapid  change  of 
T„  values  as  shown  in  Fig. (59)  is  chemical  exchange  between  the  OH 
group  of  hexanol  and  the  OD  group  of  DO  water.  The  Anderson-Fryer 
exchange,  namely  the  deuteron  exchange  between  bound  and  free  DO 
molecules,  may  be  another  possibility  but  it  would  be  hard  to  confirm. 
In  order  to  investigate  chemical  exchange,  a  sample  of  hexanol 
microemulsion  with  composition  W/0  =  0.6;  S/0  =0.2  was  prepared  with 
a  mixture  of  10%  DO  and  90%  H90  by  volume  for  high  resolution  NMR(FT) 
study.  High  resolution  proton  spectra  were  obtained  as  a  function  of 
temperature  on  a  Bruker  HX90  spectrometer.  The  spectrometer  was 
operated  at  a  resonance  frequency  of  90  »IHz.  Each  spectrum  is  the 
result  of  Fourier  transformation  of  a  large  number  of  signal  averaged 
scans. 

These  spectra  are  displayed  in  Figs . (60-64) .  The  complete  spectrum 
as  it  appears  at  80  C  is  shown  in  Fig. (60).  Only  the  protons  from  HC 
of  potassium  oleate,  the  OH  of  hexanol  and  water,  and  the  CH„  of  the 
CH„0H  group  in  hexanol  are  well  resolved  at  this  temperature.  The 
remaining  protons  contribute  the  high  intensity  structure  on  the  right. 
At  80  C  the  OH  produce  a  single  narrow  peak  and  thus  chemical  exchange- 
is  rapid  and  complete  between  the  water  and  hexanol  OH.  Figure (61)  shows 
the  relevant  left  portion  of  the  spectrum  at  23  C  on  a  250  Hz  full 
scale  sweep  [1/4  of  the  scale  of  Fig. (60)].  The  important  fact  is  that 
the  OH  of  hexanol  and  DOH  are  now  partially  resolved,  with  the  hexanol 
appearing  as  a  shoulder  on  the  stronger  DOH  peak.  The  fact  that  these 
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peaks  are  not  completely  resolved  means  that  even  at  room  temperature 
a  significant  amount  of  exchange  takes  place.  However,  the  averaging 
is  only  partial  and  this  is  therefore  a  region  of  intermediate  exchange. 
As  the  temperature  is  raised  to  60  C,  the  shoulder  disappears [Fig.  (62)  ] 
and  the  line  shape  becomes  symmetric.  This  means  that  exchange  between 
the  two  different  OH  groups  becomes  rapid  enough  to  cause  complete 
averaging  and  the  two  peaks  are  completely  overlapped.  As  the  tempera- 
ture is  further  increased  to  80  C  [Fig. (60)  and  Fig. (63)],  the  lines 
become  narrower  due    presumably  to  decreasing  viscosity  of  the  liquid. 
The  spectra  are  completely  reproducible  on  thermal  cycling.  Figure(64) 
obtained  at  23°C  only  30  minutes  after  the  80  C  spectrum  is  again 
identical  with  Fig. (61) .  These  results  clearly  demonstrate  that  chemical 
exchange  is  taking  place  and  is  changing  in  the  region  where  T  data 
show  anomalous  behavior. 

The  numerical  chemical  shift  and  line  width  data  for  the  OH  are 
given  in  Table (6).  The  shifts  are  given  in  ppm  relative  to  the  CH  OH 
peak  of  hexanol.  The  strong  temperature  dependence  of  the.  shifts  and 
of  the  linewidth  is  to  be  noted.  The  chemical  shift  decreases  by 
about  a  factor  of  four.  The  linewidth  reduction  occurs  mostly  above 
60°C. 

These  chemical  shift  data  can  be  used  to  estimate  the  lifetime 
of  the  hydroxyl  protons  bound  to  hexanol.  From  Table (6)  the  average 
chemical  shift  difference  at  23  C  between  the  two  OH  peaks  is 
approximately  0.045  ppm,  equivalent  to  4  Hz  at  90  MHz.  Therefore  the 
lifetime  of  the  OH  is  about  250  ms . 
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Fig.  (59)   The  temperature  dependence  of  spin-spin  relaxation  time 
of  the  hexanol  microenuilsions  above  room  temperature 
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Table(6)   Temperature  dependence  of  the  OH  peaks  in  the  proton 
spectrum  of  the  hexanol  microemuision  system* 


+ 
Chemical  Shifts         Full  Width  of  DOH 


Temperature 

(°C) 

+ 

2  3 

+ 

60 

+ 

80 

+ 

23 

0H(of  hexanol)   0H(of  DOH) 


peaks  at  1/2  maximum 
intensity  (Hz) 


M..25  1.215                            2.2    ±   0.3 

0.841  ^2.1 

0.654  0.8   ±   0.3 

M..27§  1.215                            2.4   ±   0.3 


*  Water/Oil  =  0.6  (  V/V  )        (  Oil  =  Hexadecane  ) 
Hexanol/Oil  =  0.4  (  V/V  ) 


3 

Potassium  Oleate/Oil  =0.2  gm/cm 


ppm  from  CH  OH  of  hexanol 


§   Shoulder  on  DOH  peak 
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IV. 2   The  Liquid  Crystals 

IV . 2 . 1   Explanation  of  Line  Spectra  and  Relaxation 

In  the  liquid  crystal  phase  a  quadrupole  splitting  was  observed 
as  discussed  in  section  II. 5  of  chapter  II.  The  quadrupole  splittings 

are  on  the  order  of  a  few  hundred  Hz,  which  is  much  less  than  the 

2 
static  e  qQ  values  of  deuterons  in  H  bonds,  showing  that  some  degree 

of  motional  averaging  is  effective.  This  motion  is  anisotropic,  i.e., 

the  D„0  molecules  maintain  a  small  degree  of  preferential  orientation, 

2 
or  else  e  qQ  would  be  zero.  The  observation  that  the  spectra  are 

powder  patterns  indicates  an  interpretation  where  the  water  molecules 
are  mostly  contained  in  lamellae   comprising  domains.  Within  these 
domains  water  molecules  may  be  "free,"  having  no  preferential  orienta- 
tion and  zero  quadrupole  interaction,  or  bound  to  the  surfactant  layer 
with  a  prefered  orientation  and  a  net  interaction.  Therefore,  the 
observed  quadrupolar  interaction   <  X  >  =  P^  X  +  P^  X^  =  P^  X^ .  The 
orientation  of  the  domains  relative  to  the  laboratory,  i.e.  the  applied 
magnetic  field,  is  random.  As  the  quadrupole  splitting  is  orientation 
dependent  a  characteristic  powder  pattern  is  observed.  The  patterns 
observed  are  characteristic  of  a  system  possessing  axial  symmetry.  The 
dimensions  of  the  domains  are  sufficiently  extensive  that  water 
molecules  evidently  do  not  diffuse  from  one  domain  to  another  during 
the  characteristic  time  of  the  NMR  experiment.  At  junctions  between 
domains  there  will  be  regions  of  disorder  where  D„0  molecules  will 
tumble  isotrcpically .  X,  is  expected  to  be  a  function  of  temperature 
decreasing  with  increasing  temperature  whereas  P   is  expected  to  be  a 
function  of  composition,  proportional  to  the  area  of  the  water /oil 
interface. 
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Figure(65)  and  Fig. (66)  show  the  wide  line  shape  of  a  typical 
lamellar  hexanol  liquid  crystal  W/0  =1.4   at  different  temperatures. 
The  two  types  of  water  have  different  relaxation  times  which  are 
temperature  dependent.  Below  the  freezing  point  of  about  4  C,  the 
formation  of  ice  crystals  evidently  disrupt  the  surfactant  layer,  and 
releases  some  of  the  trapped  water  which  then  becomes  mobile  and  gives 
rise  to  the  prominant  narrow  central  line.  Eventually  the  wide  line 
completely  disappears.  The  relatively  broad  central  line  is  actually 
a  collection  of  many  narrow  lines  which  are  overlapped  by  broadening 
at  lower  temperature.  The  central  line  will  disappear  also  when  the 
temperature  is  cooled  further.  The  unobserved  intensity  is  presumably 
spread  over  several  hundred  kHz,  characteristic  of  rigid  D„0 . 

In  a  heating  measurement  above  room  temperature  we  also  observe 
growth  of  the  central  line  due  to  thermal  agitation.  Above  a  certain 
temperature,  the  wide  line  disappears  due  to  phase  separation  and 
disruption  of  the  lamellar  structure  resulting  in  smaller  domains  with 
an  increased  probability  of  water  molecule  transfer  between  domains. 

For  further  quantitative  comparison,  we  introduce  Fig. (67)  to 

show  the  correspondence  of  the  true  line  shape  to  the  derivative  mode. 

1/2    -1 
Therefore,  we  can  calculate  T   by  using  the  formula  T  =  (tt  3    W) 

W  is  the  line  width  of  a  single  line,  Av  is  the  line  splitting  corres- 

3   2 
ponding  to  j(e   qQ)     as  shown  in  Fig. (25).  W  is  a  sum  of  intrinsic 

width  and  broadening  contributed  by  magnetic  field  inhomogeneity  W . 

From  experimental  data,  W  can  be  as  large  as  ±30  Hz.  The  calculated 

T„  varies  from  1.22  to  2.1  ms  over  a  wide  range  of  temperature  and 

this  is  comparable  with  2.1  ms  measured  by  the  spin-echo  pulsed  method 


131 


t  =  27.7  C 


■V1^ 


/       V 


\m^m^^ik\ 


328  H 


\  A** 


^^"* 


x«/ 


t  =  -9.3"C 


'Alt»*W 


o 


360  H 


t  =  -19"C 


••-— H, 


l*V 


'XyvAN*^ 


f 57  H  — H 

z 


Fig. (65)   Line  shape  variation  of  the  liquid  crystal  W/0  =  1.4  with 
S/0  =0.2  obtained  on  cooling  the  sample  from  room 
temperature  to  -19  C  (vertical  scale  is  arbitray) 
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Fig. (66)   Line  shape  variation  of  the  liquid  crystal  W/0  =  1.4  with 
S/0  =0.2  obtained  by  heating  the  sample  from  room 
temperature  to  90  C  (vertical  scale  is  arbitrary) 
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ABSORPTION  LINE  SHAPE 


DERIVATIVE  MODE 


W:  LINE  WIDTH 


Fig. (67)   The  powder-pattern  wide-line  spectra  in  the  derivative  mode 
and  the  corresponding  absorption  line  shape 
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as  shown  in  Fig. (73).  The  quadrupole  splitting  Av   is  about  200^10  Hz 
over  a  wide  range  of  temperatures .The  chief  source  of  uncertainty  in 
the  splitting  is  due  to  the  growth  of  the  central  peak  as  the  tempera- 
ture is  lowered,  which  superposes  it  on  the  wide  line  spectra. 

As  to  the  hexanol  cylindrical  liquid  crystal  W/0=1.0,  a  larger 
splitting  (about  275±10  Hz)  is  observed,  but  the  line  width  is  about 
the  same  as  for  the  sample  with  W/0=1.4  and  fluctuates  around  120-30  Hz, 
over  a  wide  range  of  temperature.  This  is  also  verified  by  the  pulsed 
method  to  have  the  same  1'   values  of  2.2  ms  [  Fig .  (71)  ]  .  The  larger  quadra- 
pole  splitting  means  less  motional  narrowing.  This  is  reasonable  because 
the  cylindrical  rods,  with  diameter  roughly  100  A,  has  more  restriction 
for  the  polar  groups  of  the  surfactants  than  does  the  lamellae  W/0=1.4. 

The  results  from  the  pulsed  method  indicate,  that  the  T„  curve  of 
the  liquid  crystals  has  thermal  hysterisis  below  the  freezing  point.  In 
the  supercooling  region,  the  liquid  crystal's  geometry  has  not  been 
destroyed  by  ice  crystal  yet,  the  short  T„  spin-spin  relaxation  time  is 
still  dominated  by  the  bound  water  contribution.  Once  the  ice  crystal 
aggregates  are  large  enough  to  disrupt  the  boundaries,  the  remaining 
water  can  travel  around  more  freely  as  mentioned  before,  and  we  observe 
a  much  larger  1,  ,  which  again  decreases  when  the  temperature  is  cooled 
further.  In  the  heating  process,  the  liquid  crystal  structure  by 
freezing  is  limited  in  local  area,  because  rapid  melting  starting  from 
a  temperature  colder  than  the  lower  limit  of  supercooling  quickly 
restores  the  same  short  T„  value.  But  this  is  not  true  in  the  case  of 
phase  separation  in  the  higher  temperature  region.  The  same  separation 
point  is  always  observed  on  heating.  But  if  we  cool  down  rapidly  from 
an  initial  temperature  above  the  phase  separation  point,  we  do  not 
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observe  a  short  T   right  below  the  separation  temperature;  instead 
it  takes  many  hours  or  a  day  for  T   to  return  to  its  original  value. 

The  liquid  crystal  W/0  =1.2  has  its  composition  between  that 
of  cylindrical  rod  and  lamellae.  This  sample  has  a  quadrupole  splitting 
Av  -  220±10  Hz  over  a  wide  temperature  range,  smaller  than  that  of 
W/0  =1.0  but  larger  than  that  of  W/0  =1.4  obtained  from  CW-measurement , 
The  line  width  of  the  main  peaks  is  about  the  same  as  in  the  other  two 
samples.  This  sample  produces  two  kinds  of  spin  echoes  corresponding 
to  two  different  T  values  in  the  temperature  range  from  -2  C  to  18  C. 
These  values  are  denoted  by  □  in  Fig. (72).  This  kind  of  two  features 
coexistance  stops  before  room  temperature,  visual  observation  provides 
no  information  on  this  region  because,  hexadecane  freezes  and  the 
sample  becomes  opaque  at  18  C.  Above  18  C,  the  spin  echo  returns  to 
one  feature. 

All  those  samples  of  the  liquid  crystal  have  significantly 
increased  T  values  after  violent  shaking  as  shown  in  Fig. (70) .  This 
is  very  similar  to  the  T„  increasing  after  freezing  because  of  local 
disruption  of  the  structure  affecting  the  motion  of  water  molecules. 
As  observed  from  CW  spectra,  shaking  makes  the  cut-off  frequency  due 
to  quadrupole  splitting  of  powder  pattern  more  apparent  as  may  be  seen 
by  comparing  Fig. (65)  with  Fig. (68)  and  Fig. (69). 

The  composition  dependence  of  the  quadrupole  splitting  and  the 
width  of  the  wide  line  are  listed  in  Table(7).  Visual  observation 
verifies  the  temperature  of  phase  separation  is  the  same  as  the 
temperature  obtained  from  the  relaxation  time  measurement. 
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Fig. (69)   The  wide  line  spectrum  of  the  hexanol  liquid  crystal 
W/0  =1.0  with  S/0  =0.2  after  shaking 
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Fig. (71)   The  temperature  dependence  of  the  relaxation  times  of 
the  hexanol  liquid  crystal  W/0  =1.0  with  S/0  =0.2 
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Fig. (72)   The  temperature  dependence  of  the  relaxation  times  of  the 
hexanol  liquid  crystal  W/O  =1.2  with  S/O  =0.2 
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Fig. (73)   The  temperature  dependence  of  the  relaxation  times  of  the 
hexanol  liquid  crystal  W/0  =1.4  with  S/0  =  0.2 
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Table (7)   The  composition  dependence  of  quadrupole  splittings  and 
widths  of  the  wide  lines 


W/0  Av  (Hz)  W  (Hz) 

1.0  275  ±  10  120  ±  30 

1.2  220  ±  10  120  +  30 

1.4  200  ±  10  120  +  30 


**  The  above  data  are  taken  at  28  C,  but  are  true  over  a  wide  range  of 
temperatures  during  which  the  liquid  crystal  structure  remains. 


IV . 2 . 2  Motional  Narrowing  of  the  Quadrupole  Splitting 

The  quadrupole  splitting  of  several  hundred  Hz  in  the  liquid 
crystal  is  much  smaller  than  the  value  estimated  by  using  the 
relational. 5. 12). 

2 

Av  =  !•(  ^JQ  )     =  2(  258.6  kHz  )  x  |=  2.42  kHz    (IV.  2.  2.1) 
4    h   obs.    4  d 

where  A  is  the  thickness  of  the  bound  water  layer  of  roughly  about 
1.25  X,  d  is  the  thickness  of  a  single  water  layer  trapped  by 
surfactant  interfaces  of  roughly  about  100  A.  This  indicates  a 
considerable  amount  of  additional  time-averaging  of  the  electric 
field  gradients  giving  rise  to  the  splittings.  Such  time  averaging 
may  be  caused  by  three  origins:  (I)  180   flip  of  the  bound  D^O 
molecules,  which  immediately  reduces  the  electric  field  gradient  to 
one  half,  and  other  modes  of  motion  such  as  vibration  etc.  , 
(II)  tumbling  of  the  water  molecules  by  thermal  mo  Lions  while  being 
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constantly  bound  and  released  from  the  binding  sites,  (III)  motion 
(within  the  relevant,  time  scale)  of  the  group  to  which  the  water  is 
bound.  The  last  possibility  is  probably  more  important  in  the  case  of 
the  lamellar  liquid  crystals  than  for  the  cylindrical  case,  because 
the  packing  of  the  polar  groups  of  the  surfactants  on  the  flat 
interface  is  not  so  tight  as  it  is  on  the  curved  interface.  The 
relative  immobility  of  polar  groups  on  the  curved  interface,  micro- 
emulsions  for  example,  has  been  demonstrated  in  a  high  resolution  NMR 
study  on  the  surfactants  by  Hansen  [41].  This  effect  explains  why  the 
cylindrical  rod  phase (W/0=1. 0)  gives  a  slightly  larger  characteristic 
splitting,  because  the  circular  geometry  of  the  cylindrical  rod  makes 
the  polar  groups  more  restricted  in  motion.  If  this  possibility 
dominates,  the  water  molecules  provide  a  useful  and  non -perturbing 
probe  for  investigating  motions  of  the  binding  groups.  If  the  second 
possibility  dominates,  the  deuteron  NMR  spectrum  should  be  determined 
primarily  by  the  energy  of  binding  of  the  water  of  hydration.  It  is 
informative  to  establish  the  set  of  relevant  time  scales  before  going 
further. 

The  criterion  for  the  start  of  motional  averaging  is  that  the 
frequency  of  the  averaging  process  be  comparable  to  the  transition 
frequency  between  the  nuclear  energy  levels.  In  a  sample  containing 

all  possible  orientations  of  the  OD  bonds,  only  a  mean  value  can 

2  n 
,  ,e  qQ, 

be  used  for  the  critical  frequency;  this  is  or  the  order  of  (    -J 
for  a  rigid  lattice  D?0  molecule,  giving  a  critical  timescale  for 
the  averaging  process  of  about  10   sec.  Averaging  processes  occurring 
over  a  longer  timescale  than  this  are  unobserved  and  one  obtains  only 
a  spatial  averaging  of  the  electric  field  gradients,  producing  the 
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characteristic  powder  spectrum.  In  order  to  determine  whether  motion 

of  the  binding  group  does  contribute  to  the  narrowing  of  the  splitting 

of  the  water  resonance,  it  is  therefore  necessary  to  consider  the 

correlation  time  of  water,  since  it  is  not  clear  which  modes  of  the 

polar  groups  are  relevant. 

If  we  take  the  liberty  to  chocse  the  activation  energy  of  the 

binding  site  as  3  kcal/mole  from  Fig. (16)  by  assuming  an  isotropic 

rotational  motion[25],  the  rate  to  disrupt  the  water  binding  group 

will  be  on  the  order  of  (~th    exp(-3  kcal/RT)  -  10   sec   .  If  any 

n 

free  water  molecule  has  its  correlation  time  shorter  than  10  "   sec, 
then  the  reorientation  between  free  and  bound  molecule  is  possible. 

The  rotational  correlation  time  of  a  free  water  molecule  is  generally 

-12 

about  10   "  sec,  thus  we  expect  the  rapid  bonding  and  releasing  from 

the  binding  sites  occurring. 

For  the  bound-free  two  state  model,  we  need  two  correlation  times 
to  characterize  the  motion.  Therefore, 


T  =  ^Tb  +  (  1  -  A  )  Tf  (IV. 2. 2. 2) 

c   d   c         d     c   , 


and  d  have  been  defined  i.e.,  A  =  thickness  of  the  bound  water  layer 
and  d  the  thickness  of  water  between  surfactant  interfaces.  Since  the 

T  values  of  the  liquid  crystal  region  in  Fig. (71-73)  are  so  different 

-2 
from  those  of  the  isotropic  case  (about  10   factor  of  the  T2  in 

microemulsions ,  which  experience  rapid  exchange  between  bound  and  free 

water) ,  it  is  reasonable  to  assume   T   >>  T 
*  c     c 

Therefore, 

T  ~  A  Tb  (IV. 2. 2. 3) 

c   d   c 
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In  the  special  case  of  rapid  reorientation  about  an  axis  which 
itself  tumbles  fairly  slowly,  only  the  tumbling  will  determine  the 
relaxation  rates.  Woessner [42] [43]  provides  the  theoretical  treatment 
in  this  special  case  and  obtains  the  equations 


(  -i—  h   ,    .   x  I  (  3  cos2'Ji  -  1  )2,         (IV. 2. 2.4) 


T   „     I,  .   isotropic   4 


where  (  T.  „  ) .   „    .   are  given  by 
1,2   isotropic 

2         2  t         8  T 

(  I  )  =  -^#) 2  ( „  C„   +  !t)   (IV .  2  . 2  .  5 ) 

(  ^isotropic   80^  <R  >    \   +  ^2  ^2  ±  +  ^2   J> 

o   c  o   c 

2  5  T         2  x 

/In  -   3,e  qQ.2       L_  + £ ).(IV.2.2.6) 

(  ^isotropic  ~  80(  h  }  (3^c  +      2  2  +       +   2   2>'<1V' 

o   c         o   c 

2  2 
Here   w  X    is  no  longer  much  smaller  than  1.  we  have  to  use 
o  c 

(IV. 2. 2.  5)  and  (IV.  2.  2.  6)  instead  of  (11.3.5).  U>   is  the  angle  between 
the  axis  of  fast  rotation  and  the  principal  axis  of  the  electric  field 
gradient  tensor  (OD  bond).  Since  the  relaxation  time  measurement  gives 
a  very  similar  result,  we  can  take  just  one  typical  case  of  W/0  =1.4 
at  28°C.  T   and  T   are  270  and  2.1  ms  respectively.  We  obtain  Tc  = 
0.32  x  10   sec  by  dividing  the  two  Woessner  equations  (IV.2.2.4)  with 

W  =  2tt  x-  6.5  x  10  radians/sec.  This  I   is  for  the  slowest  component 

Wo  _  -6 

of  the  motion  of  the  bound  water.  T   is  quite  close  to  ^  -  0.154  x  10 

second.  This  explains  roughly  why  the  quadrupole  splittings  of  the 
liquid  crystals  are  so  small. 
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I V .  3   Comparison  between  Uexanol  and  Pentanol  in  the  Water-Oil  System 

IV . 3 . 1   Electrical  Conductivity  Differences 

In  the  paper  "Some  Structural  Aspects  of  Microemulsions  and 
Cosolubilized  System"  by  Shah,  Walker,  Hsieh,  Nelander  and  Pepinsky 
(Presented  in  the  Improved  Oil  Recovery  Symposium  of  the  Society  of 
Petroleum  Engineers  of  All  IE  in  Tulsa,  Okla. ,  March  22-24,  1976)  it  is 
reported  that  the  electrical  resistance  of  microemulsions  containing 
hexanol  and  pentanol  as  the  cosurfactant  have  dramatic  differences. 
(Microemulsion  System:  hexadecane  +  potassium  oleate  +  cosurfactant 
+  water . ) 

Figure (7 5A)  and  Fig. (75B)  are  reproduced  from  this  paper.  A 
cosolubilized  solution  model  was  proposed  to  explain  the  low  resistivity 

of  the  pentanol  system.  There  is  no  explanation  for  why  the  addition 

H 
of  one  more  -C-  group  to  pentanol  causes  such  dramatically  different 

H 
results . 

Measurements  of  the  NMR  relaxation  times  T   and  T?  obtained  by 

the  spin-echo  method  provide  new  information  for  comparing  these  two 

systems.  This  evidence  is  discussed  in  the  following  subsections. 

I V . 3 . 2  Exper imental  Evidence  and  Explanation 

For  comparison,  relaxation  times  T-.  and  T   of  deuterons  were 
measured  over  a  wide  range  of  composition  in  the  oil/water  systems 
containing  hexanol  and  pentanol.  Data  at  room  temperature  are  shown 
in  Fig. (70)  and  Fig. (76).  We  observe  that  the  pentanol  system  starts 
having  liquid  crystalline  properties  (low  T  value  of  about  a  few 
miliseconds)  with  lower  water/oil  ratio  than  that  of  the  hexanol  system. 
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The  hexanol  system  is  a  perfect  microemulsion  for  the  composition 
W/0  =  0.6,  but  the  pentanol  system  is  not.  Careful  relaxation 
measurements  of  pentanol  W/0  =  0.6  were  made  from  -40  C  to  100  C. 
Below  25  C,  we  observed  an  abnormal  spin  echo  envelope.  Normally,  the 
envelope  of  the  spin  echoes  starts  at  the  peak  of  the  free  induction 
decay.  As  we  see  in  Fig. (77)  although  the  free  induction  decay  is 
large,  the  echoes  are  much  smaller  than  would  be  expected.  We  dc  not 
observe  another  feature  of  echoes  associated  with  short  T9  as  has 
been  observed  in  the  hexanol  liquid  crystal  W/0  =1.2  in  the  lower 

right  portion  of  Fig. (74) .  This  implies  the  envelope  of  the  short  T„ 

+         .  +..-1 
coincides  with  the  free  induction  decay  with  T„=  T   ,  where  (T„) 

is  approximately  the  half  width  at  half  intensity  of  the  free 

induction  decay  due  to  the  field  inhomogeneity .  The  T     value  is  about 

0.7  ras.  This  situation  with  short  T   is  similar  to  the  liquid  crystal 

behavior.  Figure(78)  shows  an  intermediate  range  between  microemulsion 

and  liquid  crystal  behavior  at  different  temperatures.  At  5  Cs  a 

strong  liquid  crystalline  tenaency  is  apparent  (a  -  a  -  1)  whereas 

near  room  temperature  liquid  crystalline  tendency  is  minor  (a  -  a  ~  0) , 

Similar  phenomena  occurs  at  lower  temperature.  In  Fig. (79)  , 

graph(a)  shows  spin  echoes  of  the  normal  type  at  -6.98  C.  Graph(b) 

and  graph (c)  show  that  the  train  of  echoes  disappear  completely  when 

cooled  down  to  -9.3  C;  the  observed  T   is  then  only  T  =  T  =  0.72  ms . 

Graph (d)  and  graph (e)  show  the  echoes  reappear  on  further  cooling 

and  reach  their  normal  type  at  -13.95  C.  These  observations  are 

reversible  on  thermal  cycling.  The  small  hump  of  spin  echoes  in  both 

graph (b)  and  graph (d)  is  believed  to  be  a  modulation  of  the  spin 

echoes  due  to  the  coupling  (dipolar  or  scalar)  between  the  resonant 
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and  nonresonant  deuteron  spins  subject  to  the  local  field  inhomogeneity 
[44]J45],[46]J[47]. 

The  results  of  visual  observation  indicated  in  Fig. (80)  are  that 
the  pentanol  dispersion  W/0  =0.6  gives  a  dark  pattern  for  all  angles 
of  the  two  polarized  plates  relative  to  each  other  at  room  temperature. 
In  the  normal  microemulsions  or  cosolubilized  molecular  solutions,  we 
expect  to  see  the  brightest  clear  pattern  when  we  orthogonalize  the 
polarized  plates  because,  the  scattered  light  has  its  biggest  intensity 
at  90°  as  we  view  it  from  the  front  side  of  the  light  source.  If  it  1S 
a  liquid  crystal  phase  with  its  grating-like  "slit"  width  comparable 
to  the  visible  light  wave  length,  we  will  see  interference  patterns 
over  a  wide  range  of  rotating  angle.  The  normal  hexanol  liquid  crystals 
with  their  "slit"  width  of  about  roughly  100  A  (  the  width  measured 
with  electronmicrograph  by  Shah  and  Pepinsky  in  the  paper  mentioned 
in  the  last  subsection  is  62  R  for  hexanol  lamellae  liquid  crystal), 
which  is  small  compared  with  visible  light  wavelength  of  several 
thousand  X,  therefore  we  can  see  interference  pattern  only  when  the 
rotating  angle  is  close  to  90°  so  that  the  path  difference  can  be 
integer  multiples  of  the  wavelength  to  give  birefringence. 

The  dark  pattern  implies  that  the  water  layer  is  so  narrow  that 
all  incoming  light  beams  can  not  penetrate  the  sample  material.  Since 
the  pentanol  dispersion  T-j/0  =0.6  has  a  very  short  1,^    (T^=  1,^) 
corresponding  to  a  rather  long  correlation  time,  which  is  proportional 
to  A/d  by  the  analysis  in  the  last  section,  the  water  layer  thickness 
of  the  liquid  crystal  phase  is  quite  small  and  is  consistant  with  the 
dark  pattern  result  by  visual  observation.  The  portion  of  the  liquid 
crystal  phase  in  the  pentanol  dispersion  W/0  =  0.6  is  very  small  so 
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that  unlike  Fig. (78)  there  is  no  difference  | a  -  a  |  ;  i.e.,  the  echoes 
are  normal.  For  temperature  above  32  C,  it  is  conjected  that  the 
thermal  expansion  and  reformation  of  the  domains  due  to  molecular 
motion  increase  the  thickness  of  the  water  layer,  only  the  light  with 
short  wavelength  is  possible  to  have  interference  pattern.  This  is 
why  we  see  very  light  purple  birefringence.  The  weak  birefringence 
disappears  again  at  72  C  and  becomes  bright  transparent,  when  the 
thermal  motion  of  the  water  molecules  are  strong  enough  to  break  the 
liquid  crystal  phase. 

Even  though  the  portion  having  liquid  crystal  characteristic  is 
not  apparent  in  the  pentanol  dispersion  W/0  =  0.6  at  room  temperature, 
it  may  still  provides  an  adequate  path  for  the  translational  motion 
of  charge  carriers  such  as  potassium  ions,  in  which  case  we  can 
understand  why  the  electric  conductivity  is  much  higher  than  in  the 
hexanol  micro emuls ions ,  which  have  isolated  water  spheres  to  prohibit 
the  motion  of  the  potassium  ions. 
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Fig.  (76)   Composition  dependence  of  the  relaxation  times  of  the 
pentanol  dispersion  system  with  S/O  =0.2 
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Fig.  (77)   The  pattern  of  spin  echo  of  the  pentanol  emulsion  W/0  =  0.6 
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Fig. (79)   The  change  of  spin  echoes  of  the  pentanol  emulsion  W/0  =0.6 
with  S/O  =  0.2  on  cooling 
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Fig. (80)   The  temperature  dependence  of  the  relaxation  times  of 
pentanol  emulsion  W/0  =0.6  with  S/0  =0.2 


157 


IV . 4   Conclusions 

Continuous  wave  and  pulse  techniques  have  both  been  used  to  study 
deuteron  NMR  lineshape,  resonance  frequencies,  and  relaxation  times 
T  and  T„  in  oil-water  dispersions  as  a  function  of  temperature  and 
composition.  A  model  in  which  the  water  molecules  rapidly  exchange 
between  a  "bound"  state  on  the  surfactant  interface  and  a  "free"  state 
in  the  bulk  liquid  is  highly  successful  in  explaining  many  of  the 
observations. 

For  the  micro emulsions  the  relaxation  rate  is  proportional  to  the 
reciprocal  of  the  radius  of  the  water  droplet  spheres  dispersed  in  the 
oil.  Analysis  of  the  data  and  calculation  of  the  relaxation  times  of 
bound  water  using  f ormula(IV. 1. 1. 6)  demonstrates  that  only  one  molecular 

layer  is  bound  to  the  polar  groups  of  the  surfactants.  The  correlation 

-12 
time  for  both  bound  and  free  molecules  is  on  the  order  of  10    sec, 

which  is  fast  compared  with  the  reciprocal  of  the  resonance  frequency 

(  10   sec  ).  The  electric  quadrupole  interaction  is  averaged  to  zero 

by  the  rapid  interchange  of  water  between  the  two  states  and  diffusion 

of  "bound"  molecules  around  the  spherical  surface  of  the  droplets. 

High  resolution  proton  spectra  obtained  on  a  partially  deuterated  sample 

(water:  90%  HO  +  10%  DO  original  mixture)  clearly  demonstrate  that 

chemical  exchange  takes  place  at  room  temperature  and  higher  between 

the  OH  of  hexanol  and  DOII.  At  2  3°C  the  hexanol  OH  appears  as  a  partially 

resolved  peak  on  the  strong  water  OK  peak  showing  therefore  a  condition 

of  intermediate  exchange.  As  the  temperature  is  raised  to  60  C  the 

exchange  becomes  complete  and  the  lines  merge  and  then  narrow  still 

further  due  to  decreasing  viscosity  as  the  temperature  is  increased  Lo 

80°C.  This  chemical  exchange  is  believed  to  be  correlated  with  an 
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anomalous  behavior  of  the  spin-spin  relaxation  time  T^   occurring  in 
the  vicinity  of  57  C. 

lhe  cylindrical  rod  and  lamellar  liquid  crystal  phases  exhibit  a 
static  quadrupolar  splitting  of  several  hundred  Hz,  which  was  most 
readily  observed  as  a  CW  powder  pattern  spectrum.  The  spectrum  is 
drastically  modified  both  at  low  temperature  due  to  disruption  of  the 
liquid  crystal  structure  by  ice  crystallization  and  at  high  temperature 
due  to  thermal  agitation  and  phase  separation.  The  observed  quadrupolar 
splitting  is  approximately  ten  times  smaller  than  is  expected  from  the 
two  state  bound-free  model  without  the  occurrence  of  additional  motional 
averaging.  The  small  value  of  the  quadrupole  splitting  is  therefore 
attributed  to  additional  degrees  of  freedom  including  (i)  librational 
and  vibrational  modes  of  the  bound  DO  molecules,  (ii)  tumbling  of 
the  water  molecules  occurring  as  bonds  with  the  polar  groups  of  the 
surfactants  are  broken  and  reformed  at  the  same  site,  and  (iii)  motion 
(within  the  relevant  time  scale)  of  the  group  to  which  the  water  is 
bound.  The  correlation  time  for  slowest  tumbling  of  the  bound  water 
at  room  temperature,  0.32  x  10~   sec,  is  only  slightly  larger  than 
the  characteristic  time  of  the  experiment,  approximately  0.15  x  10 
second.  The  quadrupolar  splitting  in  the  hexanol  cylindrical  rod 
liquid  crystal  is  slightly  larger  than  in  the  lamellar  phase.  The 
difference  is  attributed  to  the  curved  geometry  of  the  interface  which 
tightens  the  packing  of  the  polar  groups  and  therefore  restricts  their 
motion  as  well  as  that  of  the  bound  water,  thus  reducing  the  degree 
of  motional  averaging.  Vigorous  shaking  of  the  liquid  crystal  samples 
is  observed  to  substantially  increase  the  T^    This  effect  arises  from 
disruDtion  of  the  structure  which  permits  more  extensive  translational 
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motion  of  the  water.  A  similar  effect  is  observed  when  the  liquid 
crystal  long  range  order  is  disrupted  by  partial  freezing  of  the 
mixture.  Shaking  of  the  sample  results  in  a  drastic  viscosity  increase 
precisely  because  of  the  mutual  entanglement  of  the  entire  structure. 
Visual  observations  of  the  samples  as  a  function  of  temperature  show 
partial  separation  occurring  above  80  C  for  lamellar  and  60  C  for 
cylindrical  rod  structures,  in  agreement  with  the  resonance  data. 

The  pentanol  emulsion  with  W/0  =  0.6  ;  S/0  =0.2  was  surveyed 
with  relaxation  time  measurements  and  visual  observations  in  order  to 
obtain  better  understanding  of  why  it  displays  a  different  electrical 
resistivity  than  hexanol.  There  are  two  coexisting  features  of  the 
spin  echoes,  one  for  isotropic  microemulsion  and  the  other  for  a 
liquid  crystalline  structure  with  T?  =  T2  -  0.72  ms .  The  portion 
having  a  liquid  crystal  characteristic  is  temperature  dependent  and  is 
presumed  tc  provide  paths  for  charge  carriers  such  as  potassium  ions, 
which  enhances  the  electrical  conductivity.  This  is  also  supported  by 
the  results  of  visual  observation.  A  dark  pattern  is  observed  from 
18°C  to  42°C  implying  a  very  narrow  water  layer  characteristic  of  the 
liquid  crystal  phase  exists.  The  narrow  water  layer  prohibits  the 
penetration  of  the  visible  light.  Above  42  C  thermal  expansion  and 
reformation  of  the  domains  due  to  molecular  motion  increase  the 
thickness  of  the  water  layers  and  thus  gives  purple  colored  birefringence, 
which  in  turn  disappears  above  72  C  and  the  sample  becomes  clear  and 
bright.  If  the  cosolubilized  molecular  solution  were  true  for  the 
pentanol  emulsions,  the  visual  observation  should  give  the  same  bright 
and  clear  pattern  as  that  of  the  microemulsions . 
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10  REM    THTS  PRDfiRAM  CALCUI  ATFS  SPIN  LATTICE  RELAXATION  TIME  BY  AN 
13  RFM    ITERATIVE  LEAST  SQUARES  FIT  OF  THE  EQUATION  0  i  B&EXP  ( --T/T1 ) 

12  REM    TO  FIH  DATA,   AN  INITIAL  GUESS  FDR  T.1.  IS  REQUIRED. 

13  RFM    RELEVANT  REFERENCE:  J.  MAGN.  RES.  VOL.  25 »  263  (1977). 
16  REM    INPUT  DATA 

25  FILEVI-i:  "LPT:  " 

30  DIM  X$(60>  \  DIM  A (19)  \  DIM  T(19) 

3 1  D  I M  E  ( 1 9  )  \  D  I M  S  (  1 9  )  \  D  I M  f  (.19) 
35  PRINT  'SAMPLE  NAME  'i     \  INPUT  X* 
40  PRINT  "TEMPERATURE  "?  \  INPUT  X2 

45  PRINT  "NUMBER  OF  DATA  POINTS  " ?  \  INPUT  N 

5  0  P  R I N  T  c  F  I  D  A  M  P  L  I  T  U  D  E  *  »  "   T  I M  F  " 

55  POP  1=1  TO  N 

60  INPUT  A<I) 

65  PRINT  TAB (16)?"  "5  \  PRINT  PNT (27) J  *  A" » 

70  INPUT  T(I) 

75  NEXT  I 

80  PRINT  nD0  YOU  WANT  THIS  DATA  TABLE  PRINTED  ' >  VfNPUT  Y$ 

05  PRINT  "Tl  GUESS  "JMNPUT  G 

0  7  L  E  T  R  ;:-  - 1  /  G 

90  pRjNT  "DO  YOU  WANT  Tl  PRINTED  AFTER  EACH  ITERATION  SMNPU'I  Z$ 

92  PRINT*!.:  "SAMPLE  NAME  "  »X* 

9  A    P  R  I  N  T  #  1  :  n  T  E  M  P  E  R  A  T  U E  E  ~  "  J  X'  2 

96  PRINT*!.: 

98  ]T  Y$=r.»YES"  GOTO  350 

99  PRINT 

100  REM    SUMS  SET  TO  ZERO  AND  CALCULATION  BEGINS 

101  PRINT 

105  FOR  K'-O  TO  9 

110  LET  S(K)=0 

115  NEXT  K 

120  FOR  1=1  TO  N 

1 25  LET  E( I )"EXP(R*T< I ) ) 

130  LET  S(0)--S(0)+E(  I)  \  S  ( 1  )  =S  ( 1  )  +E  ( I )  "2 

135  LET  S(2)"S(2)+A( I )  \  S ( 3  )~S ( 3 )+A ( I ) *E ( I ) 

140  LET  S(4)=S(4)+A(I>*T(I)*E(I)  \  S  (  5  )  ~S  (  5  )  +T  ( I  )  *E  (  T  ) 

145  LET  S(6)"S(6)+T(  I  >*E(  I  )"2  \  S  (  7  )  -S  (  7  )  fA  (  I  )  *E  (  I  >  :KT  (  I )  "2 

150  EFT  S(S)~S(8)+E(I)*T(I)"2  \  8(9 )~S( 9 ) +E ( T ) "2*T ( I ) "2 

155  NEXT  I 

160  LET  A"--N*S(  1)  -S(0)~2  \  B  =  S  (  3  )  *S  (  0  )  -S  (  2  >  #S  (  1  ) 

!65  LET  C~S<2)*S(0)~N*S(3)  \  D-N#S ( 6 )-S ( 5 ) *S ( 0) 

167  LET  E~S ( 3  > *S ( 5 ) -S  < 2 ) *S ( 6 ) 

170  LET  fi:l.-A*S(4>+B*S(5)+C*5(6) 

175  LET  G2=U*S  (  4  )  +A:*S  ( 7  )  fE*S  (5 )  T  B*S  (  8  )  +2*C*S  (  9 ) 

ISO  LET  R=R-G1/G2 

196  IF  Z$<>"YES"  GOTO  200 

198  PRINT*!  I-.1/R 


160 


161 


200 
205 
210 

2.1.1 
215 
220 

225 
230 

235 
240 
245 
250 
2  5  5 
260 
265 
280 
290 
350 
355 
360 
365 
367 
370 
400 


IF  ABS<(Gl/62)/R><0,001  THEN  210 

GOTO  .1.0  5 

LET  M~~C/A  \  H~<S<2)-"M*S<0 )  )/N  \  Q~0 

PRINT* :l : 

FOR  1=1  TO  N 

LET  Y(I)=H+M#E(I)  \  Q-CH-  ( Y(  I  )~A  ( I )  )  "'2 

PRINT*!  :  Yd  )  .Q 

NEXT  I 

PRINT*.! 

F:'  R  J  N  T  *  1  t  "  R  E  L  A  <  A  T  1 0  N  T  I M  E  -  "  5  ■-  1.  /  E 

PR  I  NT* 1 

PRINT*! 

PR  I  NTH  I. 

PRINT*! 

PRINT*! 

PR INT* 1 

CLOSE*! 

STOP 

F;' R  I  N  T  *  .1.  J  "  F I  D  AM P L  I  T !J D E  "  f  "  T I M E  ' 

FOR  I 

PR  TNT 

NEXT  I 

PRINT*! ! 

GOTO  105 

END 


B--=  "iM 

SUM  OF  SQUARES 


FID  AMPLITUDE" 
1  TO  N 

( i  >  y  t  ( :i: ) 


APPENDIX  B 


10  PFM  HIT-  PROGRAM  OAt  CURATES  RELAXATION  TIMES  B'V 

11  REM   LEAST  SQUARES  FIT  OF  THE  EQUATION  A  +  B*EX 


AN  ITERATIVE 
P(--T/Tl)  10 

MUST  BE 

INE  300  (DO  NOT 

1  AX AT i.  ON  I [ME 


12  REM   AMPLITUDE  VS  TIME  DATA.   EXPERT MENIAL  DAD 

13  REM   ENTERED  AS  DATA  STATEMENTS  STARTING  LLiTH  I 

14  REM   EX OF ED  LINE  339).   AN  INITIAL  GUESS  FOR  RE 

15  REM   IS  REQUIRED. 

16  REM   RELEVANT  REFERENCE:  J.  MAGN,  RES,  VOL.  25v  263  (1977). 
25  FILEVtl  5  blpt:  ' 

3 0  D  I  M  X'h(60)    \    D  I M  A  (  1  9  )  \  D  I.  M  T  (  1  9  ) 

31  DIM  E ( 1 9 )  \  DIM  S ( I S )  \  DIM  Y ( 1 9 ) 

35  PRINT  "SAMPLE  NAME  "v  \  INPUT  X$ 

3 6  P R I N  r  *  I  I  "  3  A M  P  I... E  N  A M E  "  ?  X * 

37  PR  I  NT-Ill 

38  PRINT  nTin  YOU  WANT  Tl  PRINTED  AFTER  EACH  ITERATION "?  XINPUI  Z* 

39  PRINT  "DO  r'OU  WANT  DATA  TABLES  PRINTED  "»\ INPUT  Y* 

40  PRINT  "NUMBER  OF  RUNS  "f\  INPUT  P 

41  FOR  ,1-5  "!0  P 

43  PRINT  "TEMPEK'hTURE  "SXINPUT  X2 

45  PRINT  -NM-V'EP  OF  DATA  POINTS  "5  \  INPUT  N 

60  READ  A(] ') »  T(I) 

75  NEXT  I 

85  PRINT  "Ti  GUESS  "  v\ INPUT  G 

94  PRINT* J  J  'I  EMPERA'TljRF™  "  v  X2 

96  PR  INT  ir  :i  i 

98  IF  Yv>:"YFS"  GOTO  350 

9  9  PRINT 

100  REM    SUMS  SET  TO  ZERO  AND  CALCULATION  BEGINS 

1 0  i  P  R I N  T 

1.05  FOR  K--0  TO  9 

.1  1  0  LET  Si'K)=0 

115  NEXT  N 

120  FOR  1=1  TO  N 

J  25  I.  ET  E  (  I  )"::F:  XP  (R*  T  (  I  )  ) 

1 3  0  I...  E  T  S  (  0  ) :::  S  (OH  l::  (  I  )  \  S  (  .1  ) ;::  S  (  .1  )  +  F  (  I  ) '"  2 

J  35  LET  S(2)-S(2HAU:  )  \  S  (  3  )  ~S  (  3  )  4-A  (  I  )  *E  C I  ) 

j  4  0  !...  E  T  S  (  A  )  ■■■■■■■  S  (  4  )  +  A  (  I  )  *  T  (  I  )  *  E  '  I  )  \  S  (  5  )  ^ ;  '■ i  5  HTf  i  1  *  E  (  I  ) 

145  LET  S(6)=S(o)+T<  I  >*E(  T  )",2  '■-  S  (  7  »~S  (  7  )  TA  ( I  >  *L*  (  T  )  *  T  (  1  ) '"?. 

150  LET  S(8)=S(8)+E<I)*T<:n',2  X  S  (9  )=S  (9  )  +E  i  I )  "2^T<  I )  "2 

155  NEXT  I 

160  LET  A~N*S(  1  )-S(0)  "2  \  D;=S  (  3  )  *S  (  0  ) -S  (  2  )  *S  ( 1 ) 

1 6 5  L  E T  C :;" S  '■  2  )  *S  (  0  )  - N *'  S  (  3  )  \  D ■ :: N * S  (  6  )  ■- S  (  5  )  *  S  (  0  ) 

167  LET  E~S ( 3 ) *S ( 5 ) -S ( ? ) *S ( 6  ) 

170  LET  G:l  "0*8(4;  F B*S(5>  +C*S<6) 

175  LET  G2^D:t.S(4MA*S(7  >  t-E#S  (5  M  B  *S  <  8  H-2*C#S(9) 

180  LET  R=R-G1/G2 


162 


163 


196 

IF  !%<>' 

YES" 

GOTO  200 

198 

PRINT*!  J-l/R 

200 

IF  ABS< (61/G: 

!>/R><0.001  THEN  210 

201 

IF  R>0  GO TO  ; 

540 

205 

GOTO  105 

210 

LET  M=-C/A  \ 

H=<S<2)-M#S(0) >/N  \    Q=0 

21.1. 

PRINT* 1  : 

215 

FOR  1=1  TO  N 

220 

LET  Y(I)=H+Mji 

E( T  )  \  0  =  Q+(Y( I )~A( I ) )" 

2 '"'  *? 

PR I NT* 1 :y( I ) 

Q 

225 

NEXT  I 

230 

PR I NT* 1 : 

235 

PRINT*! 

"  R  E  L  <■ 

iXATION  TIME-  " f-l/R 

240 

PRINT*! 

11  A  = 

»H 

24  5 

PRINT*! 

■  B  = 

f  M 

250 

PR  I  NTH 

"SUM 

OF  SQUARE S=  " ?Q 

255 

PRINT*! 

260 

PR I NT* 1 

265 

PRINT*! 

267 

NEXT  J 

290 

CLOSE* 1 

290 

STOP 

340 

PRINT*! : "NEG< 

iTIUE  RELAXATION  TIME" 

342 

GOTO  255 

350 

PRINT*  1  :  "AhPI 

ITUDE" » '  TIME" 

355 

FOR  1=1  TO  N 

360 

PRINT*! JA(  I  ) 

-  T  (  I  ) 

365 

NEXT  I 

367 

PRINT*! : 

370 

GOTO  105 

400 

END 
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